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— '  This  paper  treats  the  problem  of  hyperbolicity ,  change  of  type  and 
nonlinear  wave  propagation  in  the  flow  of  viscoelastic  fluids.  Rate  equations 
for  fluids  with  and  without  instantaneous  elasticity  are  derived  and 
discussed.  The  equations  of  fluids  with  instantaneous  elasticity  are 
hyperbolic  in  unsteady  flow  and  can  change  type  in  steady  flow.  The  wave 
speeds  depend  on  velocities  and  stresses.  Some  estimates  of  wave  speeds  into 
states  of  rest  are  given.  For  many  of  the  popular  models  of  fluids  the 
vorticity  is  the  field  variable  which  changes  type.  The  vorticity  of  all 
fluids  with  instantaneous  elasticity  can  change  type  in  motions  which  perturb 
rigid  ones.  Experiments  and  analysis  exhibiting  vorticity  of  changing  type 
are  exhibited.  The  linearized  viscoelastic  problem  is  governed  by  equations 
having  the  properties  of  a  telegraph  equation.  The  damping  is  small  when  the 
fluid  is  very  elastic.  Elastic  fluids  have  a  long  memory,  a  large  time 
(Weissenberg  number)  for  relaxation.  The  damping  is  rapid  when  the  relaxation 
time  is  small  even  when  the  flow  is  very  supercritical.  It  is  shown  that 
steady  flow  around  a  body  is  of  "transonic"  type.  The  linearized  problem  for 
flow  over  a  flat  plate  is  reduced  to  an  integral  equation  for  the  vorticity 
distribution  on  the  plate.  The  problem  of  nonlinear  wave  propagation  is 
discussed  and  the  problems  of  nonlinear  smoothing  and  shocking  are 
considered.  It  is  shown  (^y  M.  Slemrod)  that  the  shocks  of  vorticity  can 
arise  from  smooth  data  in  some  models  and  shocks  of  velocity  in  other  models. 
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SIGNIFICANCE  AND  EXPLANATION 


Viscoelastic  liquids  exhibit  a  number  of  strange  and  unexpected  phenomena 
when  the  influence  of  elasticity  (the  "Weissenberg  number")  is  high. 

Numerical  simulations  of  high  Weissenberg  number  flows  have  encountered 
considerable  difficulties.  It  is  possible  that  some  of  these  phenomena  are 
related  to  a  change  of  type  in  the  governing  equations.  This  paper  reviews 
recent  results  relating  to  hyperbolic  structure  of  equations  for 
viscoelasticity  and  the  possibility  of  change  of  type,  which  were  obtained  by 
the  author  and  his  co-workers. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 
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HYPERBOLIC  PHENOMENA  IN  THE  FLOW  OF  VISCOELASTIC  FLUIDS 


D.  D.  Joseph 


1 .  INTRODUCTION 

The  equations  of  gas  dynamics  have  a  hyperbolic  structure  which 
support  waves  of  compression  and  rarefaction.  The  linearized  theory  of 
gas  dynamics  leads  to  the  wave  equation  in  which  compression  and  rare¬ 
faction  waves  are  on  an  equal  footing  in  that  both  types  of  waves  may 
propagate  without  change  of  form.  The  nonlinear  theory  of  gas  dynamics 
leads  to  striking  new  qualitative  understandings  of  wave  propagation. 
Initially  smooth  waves  of  compression  must  shock  up.  Initial  shocks  of 
rarefaction  must  smooth  out.  The  equations  of  steady  gas  dynamics 
change  type  from  subsonic  to  supersonic  when  the  speed  of  the  fluid  at 
some  point  exceeds  the  speed  of  sound.  The  field  of  flow  can  be  parti¬ 
tioned  into  subsonic  parts,  governed  by  an  elliptic  equation,  and  a 
supersonic  part,  governed  by  a  hyperbolic  equation.  The  problem  of 
change  of  type  in  gas  dynamics  is  called  transonic  flow. 

Quasilinear  systems  of  equations  governing  the  flow  of  incompressi¬ 
ble  viscoelastic  fluids  are  also  of  mixed  type  when  the  constitutive 
equation  admits  an  Instantaneous  elastic  response.  As  in  gas  dynamics, 
steady  flows  of  such  fluids  can  change  type. 

I  am  interested  in  the  possibility  that  many  effects  in  the  flow  of 
viscoelastic  liquids,  as  well  as  difficulties  ir  numerical  simulation  are 
associated  with  wave  propagation,  the  nonlinear  smoothing  of  shocks,  the 
shock  up  of  smooth  solutions  and  with  the  appearance  of  real  charac¬ 
teristics  and  a  change  of  type  analogous  to  the  sonic  transition  in  gas 
dynamics. 

The  notes  for  this  lecture  are  a  compendium  for  hyperbolic  things 
arising  in  the  theory  of  flow  of  viscoelastic  fluids  with  instantaneous 
elasticity. 

Fluids  with  instantaneous  elasticity  have  no  instantaneous  viscosity 
in  the  same  way  that  Maxwell  models  have  no  instantaneous  viscosity  and 
Jeffrey's  models  do.  For  example 
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is  an  Oldroyd  rate  equation  for  the  determinate  part  of  stress  t,  D/Dt 
is  an  invariant  time  derivative,  D[u]  -  i(Vu  ♦  VuT),  where  u(x,t)  is  the 
velocity,  A  is  a  relaxation  time  and  o  is  the  ',viscosity',.  Equation  (1.1) 
is  a  Maxwell  model.  It  has  a  purely  elastic  instantaneous  response.  On 
the  other  hand 

Dt  DD[u] 

*Dt  *  -  "  2n^  ♦  A-—-  (1.2) 

is  an  Oldroyd  model  with  a  retardation  time  A.  This  equation  is  of  the 
Jeffrey's  type.  The  retardation  term  produces  an  instantaneous  viscous 
response  superimposed  on  the  elastic  response.  A,  rather  than  n,  is  the 
measure  of  the  viscous  response  to  impulsive  motion.  The  most  general 
class  of  simple  fluids  with  instantaneous  elasticity  are  those  whose 
histories  are  convergent  in  the  weighted  L2  spaces  of  Coleman  and  Noll. 
Joseph,  Renardy  and  Saut  ( 1 984) ,  hereafter  called  JRS,  derived  the  gen¬ 
eral  form  of  the  rate  equation  which  is  implied  by  the  fading  memory 
theory  of  Coleman  and  Noll.  More  general  theories  of  fading  memory 
based  on  different  topologies  containing,  for  example,  models  like  (1.2), 
have  been  presented  by  Saut  and  Joseph  (1983). 

It  is  probable  that  most  polymer  solutions  have  some  instantaneous 
viscous  response  following,  say,  a  model  like  (1.2)  with  h  /  0,  but  per¬ 
haps  small.  If  the  instantaneous  viscous  response  is  small,  say  A/A  is 
small,  the  smoothing  effects  of  viscosity  will  also  be  small  leading  to 
shock  structure  of  the  type  which  is  well  known  in  gas  dynamics.  A 
theory  for  such  viscous  smoothing  of  shocks  in  viscoelastic  fluids  has 
been  discussed  by  Saut  and  Joseph  (1983)  and  Narain  and  Joseph  (1983). 
It  follows  that  the  study  of  fluids  with  instantaneous  elasticity  ought 
to  apply  to  many  real  fluids  which  have  a  small  viscous  response. 

The  visco  part  of  viscoelasticity  does  not  mean  that  the  fluid  has 
a  viscous  response.  It  means  that  the  amplitude  of  Impulsive,  shock 
solutions  will  decay  but  that  the  shocks  will  not  smooth.  This  property 
of  flow  is  not  shared  by  gas  dynamics  where  there  is  no  damping.  When 
the  damping  is  very  large  (e.g;  when  A  is  small),  the  elastic  response  is 
short.  When  the  damping  is  small  (e.g.;  when  A  is  large),  the  elastic 
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response  is  long.  In  most  problems  with  damping,  even  the  simple  ones 
which  arise  in  the  one  dimensional  theory  of  kinematic  waves  (Whitham, 

1 9711,  p.  62)  it  is  not  possible  to  form  a  shock  without  first  satisfying 
a  criterion  for  a  critical  amplitude.  Such  critical  amplitude  results 
are  now  well  known  in  the  theory  of  nonlinear  wave  propagation  in  vis¬ 
coelastic  fluids  (see  Section  13). 

There  are  possibly  many  works  on  wave  propagation  in  viscoelastic 
fluids.  Of  the  ones  known  to  me,  the  works  of  Coleman  and  Curtin  (1968) 
and  Slemrod  (1978)  seem  to  me  to  be  the  most  important.  Problems  which 
perturb  uniform  flows  of  Oldroyd  “luids  have  been  studied  by  Giesekus 
(1970),  Ultmann  and  Denn  (1970)  and  Luskin  (1984).  In  these  problems  the 
entire  field  of  flow  is  either  subcritical  or  supercritical.  There  is  no 
problem  of  transonic  type.  The  possibility  of  transonic  type  in  the  flow 
of  Oldroyd  fluids  was  first  mentioned  by  Rutkevich  (1970),  Joseph, 
Renardy  and  Saut  (1984)  treated  the  problem  of  hyperbolicity  and  change 

generally.  They  introduced  the  notion  that  in  many  models  and  in  all 

models  on  motions  perturbing  uniform  ones,  the  vorticity  is  the  hyper¬ 
bolic  variable  which  changes  type.  Too,  Ahrens  and  Joseph  (1984)  tried 
to  explain  the  striking  results  of  Metzner,  Uebler  &  Fong  (1969)  with  an 
analysis  of  the  vorticity  perturbing  irrotational  sink  flow.  This  type 
of  hyperbolic  flow  with  zero  vorticity  in  a  cone  like  region  and  nonzero 
vorticity  outside  was  observed  in  the  flow  into  a  hole. 

The  first  detailed  solution  of  flow  exhibiting  a  change  of  type  has 
recently  been  given  by  Yoo  and  Joseph  (1984).  They  consider  the  linear¬ 
ized  problem  which  arises  when  Poiseuille  flow  of  an  upper  convected 
Maxwell  model  is  perturbed  with  wavy  walls.  The  vorticity  of  this  flow 

will  change  type  when  the  velocity  in  the  center  of  the  pipe  is  larger 

than  a  critical  value  defined  by  the  propagation  of  shear  waves.  There 
is  then  a  region  around  the  pipe  axis  in  which  the  vorticity  equation  is 
hyperbolic  and  a  low  speed  region  near  the  walls  where  the  vorticity 
equation  is  elliptic.  They  linearize  the  problem  for  small  amplitude 
waviness  and  the  linearized  problem  is  solved  in  detail.  The  charac¬ 
teristic  nets  depend  on  the  viscoelastic  "Mach"  number,  which  is  the 
ratio  (M  -  U/C)  of  the  unperturbed  maximum  velocity  U  to  the  speed  of 
shear  waves  C  into  the  fluid  at  rest  and  the  elasticity  number  E.  There 
is  a  supercritical  (hyperbolic)  region  around  the  axis  of  the  pipe  when  M 
>  1.  When  M  >>  1,  the  diameter  of  this  hyperbolic  region  is  small  when 


E  is  large,  and  large  when  E  is  small.  Regions  of  positive  and  negative 
vorticity  are  swept  out  along  forward  facing  characteristics  in  the  hy¬ 
perbolic  region.  There  is  rapid  damping  of  vorticity  through  a  narrow 
layer  to  small  vorticity  in  the  core  when  M  >>  1  and  the  Weissenberg 

number  W  »  M/E  tends  to  zero.  (The  Weissenberg  number  is  proportional 
to  the  relaxation  time  of  the  fluid.) 

The  rate  of  damping  of  vorticity  decreases  as  W  is  increased. 
Flows  with  high  M  appear  to  be  more  "elastic"  when  W  is  large  in  the 
sense  that  the  damping  is  suppressed  as  the  relaxation  time  of  the  fluid 
is  increased. 

I  am  going  to  discuss  some  selected  topics  from  the  aforementioned 
papers.  I  have  introduced  3ome  new  thoughts  in  these  notes.  The  prob¬ 
lem  of  flow  around  bodies  treated  in  Section  1 1  and  the  flat  plate  prob¬ 
lem  of  Section  12  seem  not  to  have  been  considered  as  transonic  type 
problems  before.  In  Section  12  I  have  reduced  the  linearized  supercriti¬ 
cal  problem  for  flow  over  a  flat  plate  to  a  single  integral  equation.  I 
would  like  to  know  if  this  problem  can  be  solved.  Various  parts  of  the 
discussion  of  nonlinear  wave  propagation  in  Section  13  are  new. 

I  wish  to  express  my  gratitude  to  my  colleagues  and  collaborators 
in  different  aspects  of  this  work;  these  are  Michael  Renardy,  Jean 
Claude  Saut  and,  more  recently,  John  Yoo.  I  also  want  to  acknowledge 
the  help  which  I  have  received  from  my  students  Mark  Ahrens,  Edmond 
O’Donovan  and  Oliver  Riccius.  Marshal  Sleorod  has  got  a  fine  result, 
right  in  this  ballpark,  which  appears  as  Appendix  A  to  this  work.  Mark 
Ahrens  prepared  Appendix  B. 


2.  RATE  EQUATIONS  FOR  FLUIDS  WITH  INSTANTANEOUS  ELASTICITY 

It  is  perhaps  useful  to  try  to  understand  hyperbolicity  in  fluids 
with  instantaneous  elasticity  without  choosing  special  models.  We  want 
to  know  which  phenomena  are  intrinsic  and  which  are  accidental  arising 
only  for  one  or  the  other  of  the  special  models. 

A  general  theory  of  rate  equations  for  fluids  with  instantaneous 
elasticity  was  given  by  JRS.  Here  we  give  a  slightly  different  deriva¬ 
tion  leading  to  the  general  form  of  quasilinear  problems  associated  with 
motions  of  fluids  with  instantaneous  elasticity. 

T 

Let  Ct(x,t)  -  Vjl  V*  be  the  relative  Cauchy  strain  where  t) 

is  the  position  of  x.  *  (jc, t)  at  i»t.  The  determinate  stress  _t  may  &e 
expressed  by  a  functional  of  the  history  of  the  Cauchy  strain, 

t 

T  -  £  [t,Ct  (x.t)].  (2.1) 


For  the  moment  I  state  only  that  the  domain  of  F  is  equipped  with  the 
topology  of  a  Hilbert  space.  Saut  and  Joseph  (1983)  show  how  different 
choices  of  the  topological  domain  of  F  lead  to  different  types  of  con¬ 
stitutive  equations.  I  want  to  make  this  point  again,  under  (2.4),  in  the 
context  of  this  discussion  of  the  general  form  of  a  rate  equation.  F 
depends  on  the  history  of  C  (x,t)  and  explicitly  on  the  present  time, 
perhaps  through  a  kernel  function  and  in  other  ways. 

It  follows  from  differentiating  (2.1)  that 

d  t  3  -  -  dC 

dt  "  at  -  *  ]  <2-2> 

where  F,[t,£t|*]  is  a  Frechet  derivative,  evaluated  on  £t(x,t).  The  rep¬ 
resentation  theorem  for  a  linear  functional  on  a  Hilbert  space  leads  us 
to 


K(t-x,Ct(x, t)) 


<*£t 

IT  (*'T)dt 


(2.3) 
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where  K(t-i,C  (x, t))  is  a  fourth  order,  tensor  valued  kernel  function, 
not  too  singular,  satisfying  some  invariance  conditions  arising  from  iso- 

8F 

tropy  (see  JRS).  The  term  is  a  derivative  holding  the  history  of 
fixed.  This  term  is  again  a  functional  of  the  form  (2.1);  in  the  case  of 
single  integral  models  with  a  kernel  proportional  to  we  have 

3F/3t  »  -t/A  (cf*  (2.19)).  Motivated  by  this  last  observation  we  intro¬ 
duce  the  suggestive  notation 

_  def  3F 

1  ■  -‘j  ^ 

for  a  constant  A.  _t  is  stress-like  and  A  is  a  "relaxation  time". 

Saut  and  Joseph  (1983)  have  shown  that  the  smoothness  of  the  kernel 
K  in  (2.3)  depends  on  the  choice  of  the  topology  of  convergence  for  the 

dC  2 

history  C  (x,  t).  If  —  (x,*)  lies  in  the  weighted  L  (0,»)  (fading 

C  Qt  II 

memory)  spaces  of  Coleman  and  Noll,  then  the  kernel  K  must  have  at 
least  the  smoothness  which  is  required  for  the  validity  of  the  applica¬ 
nt 

tion  of  the  Schwarz  inequality  to  the  Integral  in  (2.3).  If  —  lies  in 

vl  I* 

the  (fading  memory)  Sobolev  spaces  used  by  Saut  and  Joseph  (1983),  then 

various  derivatives  of  Dirac  measures  can  appear  in  the  kernel  K.  These 

give  rise  to  "viscosity"  terms  like  the  retardation  term  in  (1.2).  The 

fading  memory  spaces  of  Coleman  and  Noll  allow  only  those  constitutive 

2 

equations  which  exhibit  instantaneous  elasticity.  The  L^(0,»)  fading 

memory  spaces  of  Coleman  and  Noll  therefore  do  not  form  a  general  basis 

for  discussion  of  problems  of  rheology.  They  do  form  a  general  basis 

for  discussion  of  problems  of  instantaneous  elasticity.  In  the  work 
-  2 

below  dom  F  *  1,^(0,®). 

Since  the  position  at  t  <  t  of  all  particles  x  on  the  same  path 
line  is  the  same 

dAt  3^t 

~dt  ’  IT  +  ^(x,t)  ‘  7lt  *  0  (2-5) 


I 


I 


I 
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we  have 

d£t  9It 

'  Vlt  +  ’  “^(u-7lt)  ♦  (u-7)7lt  -  -FtL 

where  L  -  L(x,t)  -  7u(x,t)  and  F  -  V*  (x,t)  Hence, 

™  •  t,  L>  “ 

dC.(x.t)  dF*  dF 

- -  -  F  ♦  F*  -  -L^(x,t)C  (x,  t)  -  C.  (x.T)L(x.t)  (2.6) 

at  at  — t  —  t  at  —  —  — t  —  — t  — - 

and 


<>£t 

K(t-T(Ct)^  dTJy 


dCk£ 

Kijkl  dt  dt  “  Mijkp(-,t)Lpk(-,t) 


(2.7) 


where 


W*0  •  “ 


(Kijkl  *  Kijlk)  -pi  (-'T)dT 


(2.8) 


is  symmetric  in  the  first  pair  of  indices  and  L 


(2.8)  as 


- 

pk  3x. 


.  We  may  write 


dx  i  _ 

dt  “  "  7  1 


X  - 


(2.9) 


where  M  is  the  fourth  order  tensor  valued  functional  defined  by  (2.8). 
Together  with  the  equations  of  motion 


du 

p-=  -  - Vp  ♦  div  t  ,  div  u-0, 
dt  —  ~ 


(2.10) 


8 


(2.9)  and  (2.10)  are  a  quasilinear  system  In  derivatives  of  u  and  The 
stress-like  term  ^(x.t)  and  the  coefficients  M(x,t)  are  of  lower  order 
(see  JRS)  since  they  do  not  involve  derivatives  of  u  or  and  have  other 
required  properties.  Hence,  the  system 


at 


3lij  3ue  1  - 

*  Ue  "  Mijke  3^  *  I  Tlj 


0, 


p 


3u .  3u . 

_i  ♦  u 

3t  e  3x 

e 


a 

3x/iJ 


0, 


3u 


is  linear  in  the  derivatives  of  u^,  and  p. 

'We  may  decompose  the  fourth  order  tensor  M  into  symmetric  and  skew 
symmetric  part3. 


Mijke  "  Sijke  +  \jke 

where  S  is  symmetric  and  A  is  axisymmetric  in  k  and  e.  We  note  that 

-  S.  D  , 
xjke  ek 


A.  «  . 
ljke  ek 

where 


au 

,  _ t 

’ijke  8x. 


,  i 

ljke  3x. 


ek 


2 


3u 
_ € 

ax. 


8x 


and 


f 
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Equations  governing  change  of  type  in  Oldroyd  models  arise  from 
those  Just  given  under  the  substitutions  mentioned  in  (2.19).  In  particu¬ 
lar  the  characteristic  surface  equation  becomes 


2 

a+1 

a-1  A  nl 

pu,2  - 

2 

Tn  2  t22  aJ 

♦  2[pu,u2  -  t12]<J>,<J>2  (7.12) 

♦  |pu22  -  t,,  ♦  t22  ♦  j)J  <f>22  -  0 

where  we  have  put  [Pn»P22,Alll2,A222,,Al221]  ■  [aT,,»aT22t ii2, t,2,*5( Th  -  t22)7 


8.  CONDITIONS  FOR  A  CHANGE  OF  TYPE.  PROBLEMS  OF  NUMERICAL  SIMULA¬ 
TION. 

In  general  the  condition  (7.11)  for  hyperbolicity  appears  to  be  eas¬ 
iest  to  satisfy  when  the  speed  is  large  and  the  normal  stresses  small, 
or  large  and  positive.  For  the  Oldroyd  models  (7.11)  may  be  written  as 


(1  -a 


)  -  <f.)  -  ?]  [(La)  'f  *  (a-1 ) 


♦  ( pu ,u 2  -  T,2)J  >  0 


(8.1) 


We  may  write  (8.1)  in  dimensionless  form  by  using  a  scale  length  and  a 
scale  velocity  U.  Then  dimensionless  variables  related  to  dimensional 
variables  by 


[u.v,  x.o.Y.M] 


U1  U1  T  II  111  T2j  U 

TT*  T r*  77’  77’  77’  c 


where  t#  «  r\/\  and 
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j^pu j  *  Auai  *  ^  (Pit  +  *  [2pU!U2  “  (Am*  ♦  A222l)J  gX|3x2 

+  [puj  *  A,,*,  -  1  (P„  ♦  PM)j  0  -  L.O.T.  (7.7) 

We  may  apply  the  theory  of  2nd  order  quasilinear  equations  to  (7.7). 
We  look  for  real  characteristic  surfaces  $(x,,x2)  -  const.  Let 

V$  -  (<h.$iiO).  then  (7.7)  generates  real  characteristic  surfaces  given  by 

£pui2  “  A, 22i  -  2^>>  +  +  C2pUiU2  “  (A 1 1 12  +  A2sai)]$i^t 

♦  [puj  ♦  A 1221  “  i(P„  ♦  P 2*)1  ♦»*  “  0  (7.8) 

We  note  that  if  $(x,,x2}  -  const,  gives  a  plane  characteristic  curve 

dx2  dXj  $i 

X2(x,),  then  ♦  $a  Z —  *  0.  z —  *  “  i-.  It  follows  from  (7.8)  that  we 

ax,  ax, 

have  real  characteristics  given  by 


dXj  B  ±  /Bi  -  AC 
dx,  "  A 


(7.9) 


provided  that 


B*  -  AC  i  0 


(7.10) 


where 


A  *  pu,  ■  A 1 22 1  “  2^ «»  *  ^ aa)  • 

2B  “  2pu ,u2  “  ( A 1 1 , 2  +  A222|) • 

C  »  puj1  ♦  A  I,,,  ~  2^"  +  ^a  a) 

The  condition  (7.10)  for  real  characteristics  may  be  written  as 


» 


» 
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streamlines  are  always  doubly  characteristic, 


The  streamlines  do  not  change  type,  from  hyperbolic  to  elliptic  because 
dy/dx  is  always  real  and  never  complex. 

For  problems  of  changing  type  we  have  to  analyze  the  roots  of  the 
quartic 

a  “a*  ♦  a’a,  ♦  a2a2  *aa,  ♦  a0  -  0  (7.5) 


The  roots  of  a  quartic  may  be  factored  into  the  roots  of  two  quadratics. 
The  roots  of  the  quadratics  are  real  or  complex.  Real  roots  correspond 
to  hyperboliclty  with  real  characteristic  directions. 

For  JRS  fluids  one  of  the  quadratics  has  imaginary  roots  and  the 
other  can  have  complex  or  real  roots  depending  on  the  values  of  parame¬ 
ter,  velocities  and  stresses  at  a  point.  In  this  situation  we  encounter 
a  change  of  type  with  real  characteristic  directions  only  in  certain 
regions  of  the  flow.  The  roots  which  govern  a  change  of  type  are  asso¬ 
ciated  with  the  vorticity.  The  imaginary  roots  are  associated  with  the 
relation  V2i| i  -  -5,  where  c  is  the  vorticity  and  ♦  is  the  streamfunction. 
The  streamfunction  is  associated  with  an  elliptic  operator.  Discontinu¬ 
ous  initial  data  for  the  streamfunction  will  be  smoothed. 

The  analysis  of  hyperbolicity  in  plane  steady  flows  is  easiest  when 

formed  in  terms  of  the  vorticity  C  “  Cj  •  where  c  is  governed  by  (6.4) 

in  the  plane  1,2.  After  a  little  calculation  we  find  that 

pu  u  _i!i_  .  I  P  ,Jls.  _  1  A  _J!s_  *  1  A 

J  k  3Xjdxk  2  jj  2  kjke  3Xj2xe  2  ejke  9x^xk 

-  L.O.T.,  j.k.e  -  1,2  (7.6) 


where  L.O.T.  stands  for  lower  order  terms  and  A 


ijkl 


and  skew  symmetric  in  kt.  Eq  (7.6)  may  be  written  as 


is  symmetric  in  ij 


p(au,-u2) 

0 


0 

p(au,-u2) 


-aM,m*Mlul 

_0^1112+^JJ22 

0 

au,-au2 

0 

-atM,21i  +  M,22, 

■oM12,2+M,222 

0 

0 

au,-u2 

-oM22|1  +  M222i 

-OM2212*M2222 

0 

0 

0 

a 

-1 

0 

0 

0 

(au2  -  Ua^a*®" 

♦  a*a*  ♦  a2a*  ♦ 

aia  ♦  a03  -  0 

0  au,-u2 

0  0 


where 


a*  “  Mi2l2  ~  pu  1 1 

3*  “  ~M22l2  ♦  Ml2l»  “  M»22*  +  M»U2  ♦  2  pu,u2f 

a2  *  -M22l,  ♦  M2222  ~  M12i2  -  M  l22,  +M|ui  ~  Mll22  -  p(u,+u2) 
ai  “  M222,  -  Ml2ll  ♦  Ml222  ~  M,l2(  ♦  2 pu,u2f 
a0  “  M,22i  -  pu2  . 

In  general  the  roots  of  the  6th  order  polynomial  (7.3)  are  neither  all 
real  nor  all  complex.  The  real  roots  correspond  to  real  characteristic 
lines  across  which  discontinuities  in  certain  of  the  variables  may  propa¬ 
gate. 

It  is  necessary  to  identify  the  real  characteristic  directions  and 
the  variables  to  which  they  associate.  It  is  clear  from  (7.3)  that 
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U‘1T  * 

U*1T  ' 

17  * 

m,,,217 

*  m,,2117 
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3Ti2 
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3uj 
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m,2,217 

*  "'"‘17 
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3t 

♦ 
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M22,217 
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3u, 

37 
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where  (x,,x2)  -  (x,y).  We  may  write  this  system  as 


(7.1) 


3g 

3t 


£ 


(7.2) 


where  3  stands  for  the  column  vector  whose  components  are 
[u,,U2,p.  t,i»  x2l,  x22].  The  characteristic  directions  a  =*  dy/dx  for  steady 
flows  are  determined  as  the  roots  of 


det  | Aa  -  B|  -  0. 
This  may  be  written  as 


*  V  v  V  * 


i 


I 


.-;1 
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f» 


w  II(t-T,i>n)c“,(t)dT 


p*  - 


W  T(t-T,I,n)C(T)dT 


F,(t)  -  2^  uC_I(t)  -  W  uC(t)] 

*  «!  nC(t)] 

The  last  terms  of  (6.9)  (under  the  integral)  are  of  leading  order  in 
hyperbolic  analysis  because  they  contain  first  derivatives  -  3jU^(x,t) 
of  Hl(x,t)  at  the  present  time.  These  terms  are  spoilers  because  they 
lead  to  third  derivatives  of  u,  as  in  (5.4),  which  do  not  reduce  to 
second  derivatives  of  £.  This  shows  that  not  every  popular  fluid  model 
has  a  vorticity  of  changing  type. 

It  may  be  helpful  to  draw  a  distinction  between  two  types  of  rate 
equations.  The  first  type  among  our  special  examples  depends  only  on 
the  velocity  and  the  stress.  We  may  regard  the  stress  components  as 
unknown,  to-be-determined,  field  variables,  Equations  (6.2)  and  (6.3). 
The  other  type,  like  (6.1)  and  (6.9),  has  coefficients  which  are  function¬ 
als  on  the  history  of  the  Cauchy  tensor.  One  could  say  that  these  equa¬ 
tions  are  not  "true"  rate  equations  because  the  coefficients  do  not 
depend  directly  on  the  field  variables  and  problems  based  on  these  equa¬ 
tions  must  be  posed  as  initial  history  rather  than  initial  value  prob¬ 
lems.  Rate  equations  of  the  first  type  are  clearly  convenient  because 
the  balance  equations  of  mechanics  are  also  relations  among  the  velocity 
and  stresses  and  the  whole  system  of  balance  and  constitutive  euqations 
is  then  closed.  Points  of  convenience  are  not  necessarily  points  of 
principle. 

7.  CLASSIFICATION  OF  TYPE  IN  STEADY  PLANE  FLOW 

In  plane  flows  there  are  six  equations  for  six  unknowns;  u,,  u2, 

Th,  p.  The  quaslllnear  equations  which  governs  the  flow  of 

fluids  with  Instantaneous  elasticity  may  be  written  as 


1 
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T  -  {-W  TT(t-t,I,II)  C_,( t)  ♦  W  T(t-T,I,II)C(T))dT 


where  W  -  3W/3I,  etc.  The  determinate  stress  vanishes  on  the  rest 
state  (OC-,-1);  hence 


Wn(t-T.3,3)  -  W(I(t-T,3,3).  (6.7) 

The  rate  equation  satisfied  by  2  nay  6e  obtained  from  (6.4)  by  dif¬ 
ferentiating  with  respect  to  t,  using 


dC  t 

-=  -  -L1(t)C(  t)  -  C(  T)L(t), 
at  —  —  —  — 


-  L(t)C-,(  t)  +  C-l(  t)LT(t) 
at  —  —  —  — 


S  “  tr  -  "2  tr[C( r)L(t)] 


HTT  dC”*  T 

~  -  tr-~-  -  2  tr[C-*(T)LT(t)] 
dt  dt  —  — 


(6.8) 


We  find  that 


di  1  -  t 

dl  *  i  i  -  k  £»  *  hk  * 


♦  i£>  +  {F,( x)tr[C( x)L(t)] 


*  £*(  i)tr[C_,(  t)LT(t)]}d  t 


« •  [‘.ft*  s-  ■  fe]- 


where 


»[v«  *  2*'7Vk  *  Vj'k.eJ  *  5(,"')'kW!l,rl  i]q.. 

~2(a4'),npck,mp”Ti7"tk  "  "TckeJTJm,me 

*  *[<H  *  u-»)<i-'uk)  -  (9tu-V)Ck  -  Vj.e'kJ 

*  EkeJ  \_~T~  Tmp,a  sk,p  *  ~  ’jq.e  lcurl  S>q] 


r  O  ’  l  j  t 

*  ekejL  2  Tmp,euj,mp  *  Tmp,meuJ,p  TJm,mpLlp,e; 


+  ~2~  (Tmp,eup,jm  *  Tmp,meup,j  +  TJp,meup,m 


♦  t  u  -*■  -r .  u  ♦t.  u  )] 
jp,m  p,me  Jp.me  m,p  Jp,m  m,pe 


Moreover 


'kej’jm.me  ’  lourl  dly  i>k  '  »[ourl  3?]k  ‘  *  <E'7)uk]  <6’5) 

Potential  flow  with  c«0  can  be  a  solution  of  (6.4)  only  if  the  last 
bracket  in  (6.4)  vanishes  u-7*. 

The  BKZ  model  is  a  single  integral  model  which  is  motivated  by  the 
nonlinear  theory  of  elasticity.  Let  C(t)  and  C“‘(t)  be  the  right  relative 
Cauchy  strain  tensor  with  invariants 

(I, II)  -  (trC( t),trC-1( t}) 

and  "strain  energy" 


W(t-T,I.II). 


The  determinate  stress  is  given  by 


The  upper,  lower  and  corotational  Maxwell  model  correspond  to  a-1,-1,0, 
respectively. 

We  may  say  that  the  JRS  models  are  of  the  Giesekus  type  when 


t  -  t  ♦  at 


with  mobility  constant  a.  The  model  actually  proposed  by  Giesekus 
(1982)  has  retardation  time  on  the  right  side  which  we  put  to  zero  for 
instantaneous  elasticity,  at*  is  added  to  the  left  side  and  a  -  1 


T  d  t  t] 

U—  -  Vu  jt  -  _t  Vu  ♦  x,  ♦  or2  -  2nD 


(6.3) 


Equations  of  the  form  (6.1),  (6.2)  and  (6.3)  permit  hyperbolic  waves 
of  vorticity.  The  last  term  on  the  left  of  (5.^)  may  be  expressed  as  a 
second  order  derivative  of  the  vorticity 


ekeJMJmpquq,pme  ’  2ekeJPjq  ^curl  -^q.e 

,  ~P  r  +  —  r  e  A  c 

2  mp^k.mp  2  kej  pqr  jmpq  r,me 

Then  (5.1*)  may  be  regarded  as  a  second  order  equation  for  the  vorticity. 

The  vorticity  equation  for  the  three  parameter  family  of  Oldroyd 
models  can  be  written  as  follows. 
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p[atteVc  *  *  ueujck,ej]_tklJMJmpquq,pmt 

.p[(3t  •  H*V>(i-V>uk  -  <atu-?>tk  -  ueuJ  e5kJ 

♦  c  fM  n  ♦  (M  u  ) 

kej  Jmpq.e  q,po  kjmpq,m  q,p  ,e 

~Up.eTJq,pq  4  Cup.q\jq,p\e  +  XTJm,me^ 


(5.4) 


j  _  ^  3 

_  ,  .  r —  (•)•  Equation  (5.4)  reduces  to  a  second  order  equa- 

where  (*),m  -  ^  H 

m 

tion  for  the  vorticity  plus  lower  order  terms  whenever  the  last  term  on 
the  left  of  (5.4)  can  be  expressed  In  terms  of  second  order  derivatives 
of  £  (see  (6.3)).  Potential  flow  is  possible  only  when  (5.4)  is  satisfied 
with  c-0,  u-V*. 


6.  SPECIAL  MODELS 

Special  constitutive  models  of  fluids  with  instantaneous  elasticity 
arise  from  (2.9)  by  choosing  special  forms  for  the  fourth  order  tensor  M 
and  the  stress-like  term  t.  In  fact  £  -  t  for  all  single  integral 
models  with  separable  kernels  in  exponential  form  with  time  constant  X. 
We  will  consider  four  special  families  of  constitutive  models. 

(1)  JRS  models  (introduced  by  JRS  (1984)). 

(2)  Three  parameter  family  of  "Maxwell"  models  introduced  by 
Oldroyd. 

(3)  Models  of  the  Giesekus  type. 

(4)  BKZ  models. 

The  JRS  models  arise  when  the  fourth-order  tensor  M  is  expressed  by  the 
most  general  expression  which  involves  second  order  tensors.  This  leads 
to 


djt 

dt 


PD  ♦  DPT  ♦  1(AQ  ♦  (All)7)  ♦  jx 


(6.1) 


for  an  arbitrary  second  order  tensor  P;  the  fourth  order  tensor  A  has 

*”  ~  T 

already  been  Introduced.  Coleman's  thermodynamics  implies  that  P  •  P  . 

The  Oldroyd  models  with  three  constants  arises  from  (6.1)  when 
The  Oldroyd  models  with  three  constants  arises  from  (6.1)  when 

£  ■  *1  * -1  S  a  S  1,  t  ■  i  and  A  chosen  so  that 


dient  u  -  V*.  Then  £  -  curl  u  -  0.  The  velocity  u  determines  the  path 
lines  (x,  t)  and  the  past  values  of  the  relative  Cauchy  strain  Ct(x,  t). 
This  gives  rise,  through  (2.11),  to  the  values  t(x,t)  of  the  extra  stress 
at  the  present  time  t.  It  is  not  guaranteed  that  this  computation  will 
lead  to  a  T(x,t)  such  that  curl  div  2  -  0  at  each  and  every  point  at 
which  i  “  0. 

Whenever  u  -  V*  there  are  two  more  equations  than  unknowns.  Evi¬ 
dently  only  certain  special  potential  flows  are  compatible  with  (5.1). 
Uniform  motions,  say  u  -  exU,  for  constant  U,  have  2*0  and  are  always 
potential  flows. 

Another  form-  of  the  vorticity  equation  is  fundamental  in  the  study 
of  hyperbolicity  and  change  of  type.  We  first  apply  the  substantial  de¬ 
rivative  ~  to  (5.1).  Then  we  apply  ^  curl  to  (2.10),  eliminate  curl 
dT 

div  —  and  find  that 
dt 

d*i 

p^pr  -  curl  div  (ML?)  -  0,  +  02  +  0,  (5.2) 

where 

L  -  VU.  (Vufy  -  ijUj 
£  -  curl  u, 

©i  *  -  j  curl  div  2, 


—  (curl  div  x) 
dt  — 


dt 

curl  div  -fr, 
dt 


(5.3) 


d  du  d 

'  pdt  (curldt  "  dt  Curl 


Third  derivatives  of  u  are  all  on  the  left  side  of  (5.2).  The  right  side 
contains  second  derivatives  of  u  and  2i  at  most.  (02  is  essentially  a 
second  derivative  of  t.)  Equation  (2.9)  shows  that  2  and  u  are  of  the 
same  differential  order,  so  that  the  left  side  has  third  order  deriva¬ 
tives  of  u  and  the  left  side  has  second  order  derivatives  of  u,  at  most. 


in  50 %  glycerin  and  49  S  water  solution)  using  a  rotating  cylinder  appa¬ 
ratus.  For  Maxwell  models  with  a  single  relaxation  time  we  have 
G(0)  -  p/A,  where  n  is  the  viscosity  and  A  is  the  relaxation  time.  There 
are  some  papers  which  report  viscosities  and  relaxation  times  which  are 
obtained  by  fitting  data  with  some  impusle  experiments  (for  example,  see 
Fapanastasiou,  Scriven  and  Macosko  0983),  Dodson,  Townsend  and  Walters 
(1971),  Laun  (1978).  This  type  of  rheometry  suffers  from  two  defects. 
Usually  a  special  form  for  the  kernel  in  the  linear  theory  is  assumed. 
The  kernel  is  represented  by  an  exponential  or  by  a  sum  of  exponentials. 
A  more  serious  defect  is  that  it  is  usually  assumed,  wrongly,  that  the 
early  time  response  to  impulse  which  gives  G(0)  satisfies  an  overly  sim¬ 
plified  asymptotic  (for  large  time)  theory  (see  Narain  and  Joseph  0  984) 
for  a  full  discussion  and  remedy  for  this  defect).  From  the  impulse 
data  reported  in  many  different  papers  we  have  estimated  that  wave 
speeds  for  water-based  polymer  solutions  are  of  the  order  of  10  cm/sec. 
The  wave  speeds  of  polymer  melts  are  of  the  order  of  100  cm/sec.  At 
the  same  time,  the  large  viscosities  of  melts  reduce  the  possibility  of 
achieving  large  velocities.  Hence  for  the  flow  of  melts  we  expect  that 
the  flow  is  usually  subcritical  M  <  1. 

We  (I,  Oliver  Riccius  and  Gordon  Beavers)  have  developed  a  wave 
speed  meter  which  appears  to  give  direct,  reproducible  values  for  C.  A 
patent  has  been  applied  for. 

5.  VORTICITY 

Some  important  special  models  give  rise  to  a  hyperbolic  vorticity 
equation.  It  is  of  interest  to  derive  the  equation  for  the  vorticity 

without  making  assumptions  other  than  that  the  fluid  has  instantaneous 
elasticity,  i.e.,  the  quasilinear  system  (2.11)  governs. 

The  vorticity  £  satisfies 

3 

p(^£  £  +  U’ V_j)  -  p£*Vu  ♦  curl  div  j.  (5.1) 

Equation  (5.1)  is  satisfied  by  £  -  0  if  and  only  if  curl  div  _r  -  0.  This 

condition  is  always  satisfied  by  Newtonian  fluids  for  which 

T 

x  *  n(Vu  *  Vu  )  for  a  constant  viscosity  n  giving  rise  to 
curl  dlv^  ■  nV2£.  In  this  case  potential  flow  is  always  a  solution  of 
(5.1).  More  generally  we  may  verify  that  potential  flow  u  -  V$, 

V  ♦  C'(x)  is  a  solution  of  (5.1)  if  and  only  if  div  1  is  the  gradient  of 

a  potential.  Given  the  potential  ♦  we  may  obtain  the  velocity  as  a  gra- 
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There  may  of  course  be  more  than  one  wave  speed.  If  C  is  not  real  there 
will  be  rapidly  growing  solutions.  The  reader  may  see  more  details  of 
this  type  of  analysis  in  the  studies  by  Rutkevich  (1969,  1972)  of  Oldroyd 
models. 

4.  WAVE  SPEEDS  II,  PHYSICAL 

The  subject  of  gas  dynamics  would  be  in  deep  trouble  if  it  were  not 
possible  to  know  the  speed  of  sound.  This  type  of  deep  trouble  3hould 
be  a  concern  in  the  study  of  hyperbolicity  in  the  flow  of  viscoelastic 
fluids.  In  fact  the  way  to  know  if  a  fluid  has  instantaneous  elasticity 
is  to  show  it  has  a  wave  speed.  Now  showing  that  a  certain  fluid  has  a 
wave  speed  is  not  easy.  First  of  all,  there  are  many  fluids  for  which 
there  are  no  wave  speeds.  These  common  cases  are  for  fluids  like 
Newtonian  one3,  or  even  viscoelastic  ones  with  a  big  viscous  response 
(3ay,  A/ A  is  not  small),  which  do  not  have  an  Instantaneous  elastic  re¬ 
sponse,  or  a  predominantly  elastic  response  to  impulsive  data.  Secondly, 
the  fluids  which  admit  a  wave  speed  admit  infinitely  many  wave  speeds. 
Just  as  the  speed  of  sound  in  a  gas  at  a  point  depends  on  the  thermody¬ 
namic  condition  there,  a  wave  in  a  viscoelastic  fluid  at  a  point  depends 
on  the  velocity  and  the  state  of  stress  there.  So  each  problem  gives 
rise  to  its  own  field  of  wave  speeds. 

The  simplest  type  of  wave  propagation  is  into  a  region  at  rest. 
The  wave  speed  for  this  is 

where  G(0)  i3  the  instantaneous  value  of  the  shear  relaxation  modulus 
G(s).  In  fact  all  kind3  of  small  amplitude  discontinuities  will  propa¬ 
gate  with  this  speed.  Jumps  in  acceleration  (Coleman  and  Curtin,  1968), 
jumps  in  velocity  and  even  in  displacement  (Narain  and  Joseph,  1982). 
The  same  propagation  speed  (4.1)  holds  for  waves  propagating  into  regions 
undergoing  rigid  motions.  We  shall  show  that  flows  which  perturb  rigid 
motions  go  supercritical  when  U/C  >  1.  This  wave  speed  is  a  fundamental 
material  parameter.  There  are  at  present  no  rheometers  for  measuring 
G(0)  (or  C)  and  there  are  no  tables  of  values.  E.  H.  Lieb,  (1975),  in  his 
Ph.D.  thesis,  measured  C  -  8.0  cm/sec  in  one  fluid  (carbomethylcellulose 
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can  be  expressed  in  terms  of  the  vorticity  £  -  curl  u.  Bernard  Coleman 
has  shown  that  his  thermodynamics  imply  that  the  fourth  order  tensor  S 
is  derivable  from  an  energy  and,  as  a  consequence 


Sijke  ”  Skeij 

3.  HAVE  SPEEDS  I,  THEORETICAL 

The  qua3ilinear  system  is  called  evolutionary  if  the  Cauchy  initial 
value  problem  for  perturbation  of  arbitrary  motion  is  well  posed  (see 
Gelfand,  1963).  We  look  at  periodic  initial  data  whose  period  is  so 
small  that  _r,u,p  are  essentially  constant  over  the  period  2ir/k,  k-|k| 
over  which 


(S^.du.flp)  *  (c.v.ir)  exp  ik*x 

are  rapidly  varying,  k-*«  .  We  seek  propagating  waves  or  the  form 

(fi_r,6u,5p)  -  (jo.v.x)  exp  i(k*x*wt)  (3.1) 

and  (jff.v.ir)  depend  only  on  the  components  of  u  and  M  which  appear  as 
coefficients  of  derivatives  in  the  quasilinear  system  (2.11).  The  result¬ 
ing  equations  are  then  divided  by  k  and  the  right  side  of  (2.11)  vanishes 
as  k  ■»  *.  This  leads  to  the  homogeneous  equations  for  the  components  of 
def 

(o,v, ir)  with  k/k  -  n  fixed 

p[“  *  h'sK  *  v  -  Vu  ■  0(i?’ 

[k  *  -  a]0lj  '  Mijkevkne  '  0(I> 

v  n  -  0 
e  e 

The  homogeneous  system  corresponding  to  (3.2)  can  be  solved  if  and  only 
if  the  determinant  of  the  coefficient  vanishes.  The  system  is  evolution¬ 
ary  if  there  are  real  wave  speeds 


Is  the  speed  of  a  wave  of  vorticity  into  a  fluid  in  uniform  motion.  The 
criterion  (8.1)  for  hyperbolicity 


[mV  ♦  |(1-a)  -  |(1*a)  -  ij  [(1  *a)|  ♦  (a-1)|  ♦  1  -  mvJ 

♦  [M*uv  -  t?  >  0  (8.3) 

depends  only  on  the  viscoelastic  Mach  number  M  (and  the  parameter 
ae[-1,1]).  When  a  -  1  (upper  convected  Maxwell  model),  the  criterion 
(8.3)  can  be  written  as 

M^v*(1*o)  +  u*(1  *y)  -  2  tuv] 

-  (1+Y)(1*a)  ♦  ta  >  0  (8.U) 

This  criterion  is  satisfied  for  high  speed  flow  if  the  stresses  o.Y.t  are 
not  large  and  of  the  wrong  sign. 

The  vorticity  equation  changes  type  at  points  at  which  the  right 
side  of  (8.1)  or  (8.3)  vanishes.  Numerical  simulations  of  problems  of 
changing  type  can  be  very  difficult.  For  example  numerical  simulations 
of  transonic  flows  have  only  become  satisfactory  in  recent  years  and 
they  are  still  being  improved.  Problems  of  the  transonic  type  can  be 
expected  in  numerical  simulations  of  the  flow  of  viscoelastic  fluid. 

Steady  flow  of  the  Oldroyd  fluids  (6.2)  is  governed  by  the  follow¬ 
ing  system  of  dimensionless  equations 

div  u  «  0, 

M*(u*7)  u  ♦  Vp  -  div  _r  -  0, 

u-Vjt  -  a(Dx  ♦  tD)  ♦  tQ  -  Or  -  D  «  -_t/W  (8.5) 

where 


W  -  UA/d 


Is  the  Welssenberg  number,  which  may  be  regarded  as  a  dimensionless 
relaxation  time  measured  in  units  of  d/U.  Stresses  relax  slowly  when  W 
is  large.  The  fluid  has  a  long  memory. 

The  high  Welssenberg  number  problem  is  that  people  doing  numerical 
analysis  can't  get  answers  when  the  stress  levels  are  high.  This  prob¬ 
lem  is  evidently  independent  of  the  "transonic"  problem  of  change  of 
type.  Most  of  the  computations  in  which  one  encounters  this  problem  are 
for  flows  without  inertia,  p  -  0.  Yoo  and  Joseph  (198*1)  also  encountered 
this  problem  in  high  speed  flows  involving  change  of  type.  They  found 
that  regions  of  positive  and  negative  vorticity  decay  rapidly  when  W  is 
small  (even  if  the  "Mach"  number  is  large).  When  W  is  large  the  super¬ 
critical  regions  are  more  "purely"  hyperbolic  in  that  the  damping  of  the 
vorticity  is  suppressed.  We  have  something  like  a  telegraph  equation 
with  damping  proportional  to  1 /W.  (See  (12.19). 

Though  the  high  Welssenberg  number  problem  and  problems  of  change 
of  type  are  Independent  they  may  be  related  problems  in  certain  special 
cases.  Though  flows  without  inertia  will  not  ordinarily  change  type,  the 
criterion  for  change  of  type  wherein  B2  -  AC  changes  sign  can  be  inter¬ 
esting  even  when  p  -  0. 

For  the  upper  convected  Maxwell  model  (a  -  1)  we  find  from  (8.1) 
that  when  p  -  0 


B2  -  AC  -  t,2,  -  [t„  ♦  j][t„  *  j]. 


The  equations  for  steady  flow  are  elliptic  when  this  is  negative.  In 
fact,  this  quantity  appears  to  approach  zero  very  rapidly  as  X  (or  W)  is 
Increased  in  numerical  integrations  of  steady  flows.  In  fact,  B2  -  AC  is 
always  negative  for  the  upper  convected  Maxwell  model 


i-  fi 


-U-t)/A  „-i. 
e  C  (  t  )d  t . 


Since  C  1  has  positive  eigenvalues,  the  principal  values  t , ,  ♦  *  T 

A  A 

are  positive  and  therefore 
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In  numerical  Integrations  using  finite  elements  or  finite  differences, 
values  of  ty  and  are  obtained  from  discrete  steps  In  which  the  cri¬ 
terion  (8.8)  may  be  violated.  In  doing  numerical  Integrations  one  should 
verify  that  the  condition  (8.8)  is  not  violated. 

Some  constitutive  models,  like  the  corotatlonal  one  with  a  -  0  may 
undergo  real  change  of  type  even  when  p  -  0. 

People  doing  flow  computations  for  viscoelastic  fluids  are  able  to 
go  to  higher  Weissenberg  numbers  when  they  have  constitutive  equations 
with  more  Newtonian  viscosity  (non-zero  retardation  times).  This  obser¬ 
vation  suggests  that  viscosity  methods  for  dealing  with  problem  of 
change  of  type  and  shocks  could  also  be  useful  for  solving  the  high 
Weissenberg  number  problem.  So  far  however  the  people  doing  numerical 
works  have  not  used  the  viscosity  method  in  the  limits  of  small  viscos¬ 
ity. 

9.  LINEARIZED  PROBLEMS  OF  CHANGE  OP  TYPE 

Up  to  now  our  3tudy  has  been  exact  and  fully  nonlinear.  We  may 
advance  our  understanding  of  the  problems  of  vorticity  of  changing  type 
by  considering  simpler  problems  which  arise  under  linearization. 

We  first  identify  a  class  of  motions,  say  motions  in  the  plane  or 
axisymmetric  motions.  We  find  some  special  exact  solution  of  all  the 
equations  which  fall  in  the  given  class.  Usually  this  special  solution 
is  featureless,  like  uniform  flow,  flow  into  a  sink,  flow  in  a  channel, 
flow  between  cylinders,  extensional  flow  and  so  on.  The  featureless 
solution  does  not  exhibit  the  unusual  features  of  change  of  type.  We 
then  perturb  all  the  equations  around  the  special  one  with  perturbations 
in  the  given  class,  and  we  linearize.  This  leads  us  to  linearized  prob¬ 
lems  with  variable  coefficients  depending  alone  on  the  special  solution 
and  not  on  the  perturbation.  We  analyze  the  linear  perturbed  problem 
for  hyperbolicity  and  change  of  type.  The  equations  for  the  characteris¬ 
tics  are  given  by  (7.10)  where  A,B,C  are  evaluated  on  the  unperturbed 
special  solution. 

The  procedures  of  linearization  are  such  that  the  characteristic 
directions  are  the  same  for  each  and  every  linear  problem  perturbing  the 
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special  one.  The  characteristic  surface?  for  linearized  problems  are  a 
gift,  since  nothing  beyond  the  special  solution  is  needed  to  compute 
them.  We  can  know  which  are  the  elliptic  regions  and  which  the  hyper¬ 
bolic  regions  of  flow  cheaply,  but  to  find  other  properties  of  the  per¬ 
turbed  flow,  say  isovorticity  and  streamlines,  we  are  obliged  to  solve 
linear  PDE's. 

JRS  (1981!)  treated  the  linearized  problem  for  a  change  of  type  in 
shear  flow  and  extensional  flow  for  the  family  of  Oldroyd  models  char¬ 
acterized  by  the  parameter  a  [1,1].  They  also  treated  sink  flow  in 
the  plane  and  circular  Couette  flow  of  an  upper  convected  Maxwell  model 
(a-1).  They  identified  the  regions  of  subcritical  (elliptic)  and  super¬ 
critical  (hyperbolic)  flow,  but  they  did  not  compute  the  characteristics 
or  solve  some  boundary  value  problem.  Shear  flow  u  -  tcy,  tc  is  the  rate 
of  shear,  is  hyperbolic  outside  a  strip  centered  on  y  -  0.  Extensional 
flow  (u,v)  -  s(x,-y),  s  is  the  rate  of  extension,  is  either  hyperbolic 
outside  an  ellipse  (where  u  and  v  are  large)  or  inside  a  region  bounded 
by  branches  of  a  hyperbola  (x2/A2)  -  (y2/B2)  <  1  (where  the  velocity  need 
not  be  large.  Sink  flow  with  radial  velocity  -  -Q/r  (potential  flow) 
is  hyperbolic  when  the  radius  r  is  small 


provided  that  the  source  strength  Q  >  2n/p,  where  n  i3  the  viscosity. 
Couette  flow  outside  a  rotating  cylinder  of  radius  a  is  hyperbolic  in  an 
annulu3  whose  inner  radius  is  either  a  or  i3  greater  than  a,  depending  on 
conditions.  When  the  inner  radius  of  the  annulus  of  hyperbolicity  is 
greater  than  a  we  have  another  example  where  the  region  of  high  speed 
flow  i3  elliptic  and  regions  of  lower  speed  are  hyperbolic.  The  varia¬ 
tion  of  the  stresses  is  important. 

Yoo,  Ahrens  and  Joseph  ( 1 98^1 )  have  tried  to  explain  some  striking 
experimental  results  of  Metzner,  Uebler  and  Fong  (1969)  with  an  analysis 
of  the  vorticity  perturbing  irrotatlonal  sink  flow.  The  main  point  at 
issue  in  the  experiments  of  Metzner,  et  al.,  is  that  they  observe  a  coni¬ 
cal  region  of  zero  vorticity.  Outside  this  region  the  flow  is  rota¬ 
tional.  If  we  accept  the  experimental  results  at  face  value,  we  must 
conclude  that  there  are  surfaces  across  which  some  derivatives  of  the 
vorticity  are  discontinuous.  This  type  of  behavior  says  "look  for  hyper¬ 
bolicity  and  change  of  type."  Such  a  field  could  not  be  supported  by  an 


elliptic  vorticity  field. 

In  the  experiments  of  Metzner,  Uebler  &  Fong  0  969)  a  fluid  is 
sucked  from  a  pipe  of  large  diameter  through  a  sudden  contraction.  If 
the  hole  into  which  the  flow  goes  is  small  the  problem  may  be  thought 
to  be  a  form  of  sink  flow.  Because  there  are  boundary  walls,  the  flow 
through  a  sudden  contraction  is  not  a  sink  flow  in  a  strict  sense.  We 
shall  imagine  first  that  the  flow  into  the  hole  is  not  strongly  influ¬ 
enced  by  the  walls  of  the  large  pipe.  We  then  have  a  hole  in  the  semi 
infinite  region  above  a  plane.  This  flow  is  then  regarded  as  an  axisym- 
metric  perturbation  of  sink  flow  without  boundaries.  The  characteristic 
surfaces  for  the  vorticity  of  all  axi3ymmetric  linearized  problems  per¬ 
turbing  sink  flow  can  then  be  obtained  by  integrating  the  differential 
equations  for  the  characteristics.  Yoo,  Ahrens  and  Joseph  then  tried  to 
determine  if  the  characteristic  surfaces  computed  in  this  way  could  be 
the  locus  for  the  discontinuity  in  the  derivatives  of  the  vorticity 

observed  in  the  experiments.  It  is  Important  to  verify  that  the  charac¬ 
teristics  are  cone  like  in  the  region  where  potential  sink  flow  was 

observed  and  that  the  region  of  potential  flow  is  in  the  spherical 

annulus  of  hyperbolicity.  These  issues  are  addressed  below. 

Metzner,  Uebler  and  Fong  (1969)  consider  high  speed  flow  of  viscoe¬ 
lastic  fluids  into  a  sudden  contraction.  They  say  that  "A  tentative 

analysis  of  the  observed  velocity  field  suggests  the  flow  upstream  of 
the  small  duct  to  be  radially  directed  toward  the  origin  of  the  spheri¬ 
cal  coordinate  system.  If  this  is  so  the  continuity  equation  gives 

ur*  -  f ( 6) .  (1)" 

They  actually  measure  velocities  in  the  cone  and  they  report  that  their 
measurements  were  accurate  and  that  f(0)  may  be  taken  as  constant  when 
0165  10#.  They  also  write  that  uQ  -  Ug  -  0  in  the  cone.  Outside  of 
the  cone  there  i3  secondary  motion  and  nonzero  vorticity. 

The  nature  of  the  comparison  of  theory  and  experiment,  explained  in 
the  captions  of  Figs.  9.1  and  9.2,  is  discussed  in  greater  detail  in  the 
paper  of  Yoo,  Ahrens  and  Joseph  (1984).  It  turns  out  that  the  regions  of 
hyperbolicity  are  model  sensitive.  Using  the  measured  values  of  physical 
parameters  in  the  equation  one  finds  that  the  observed  potential  flow 
does  lie  in  the  hyperbolic  region  (some  points  are  outside  the  hyperbolic 
region  in  Fig.  9.2(b)).  The  characteristic  surfaces  are  cone  like  in  the 
regions  where  potential  flow  was  observed. 
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Metzner,  Uebler  &  Fong  (1969)  measured  potential 
flow  in  the  sectorial  box  designated  S.  The  cross- 
sections  of  characteristic  surfaces  of  revolution 
which  are  tangent  to  the  cone  of  semi  vertex  angle 
10°  at  the  origin  are  called  C.  There  are  two  such 
surfaces  of  revolution. 
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Fig.  9.2.  Schematic  diagram  of  flov  into  a  sudden  contraction 
(Metzner,  Uebler  and  Tong  (1969)).  The  vortlcity 
appears  to  vanish  In  a  cone  with  nonzero  vortlcity 
outside.  Accepting  this,  a  jump  in  some  derivative 
of  the  vortlcity  on  the  cone  is  required.  The 
measurements  verifying  potential  flov  were  taken  in  a 
cone  with  6  *  10*  and  at  a  certain  value  in  and 
out  of  the  pipe.  The  regions  where  potential  flov 
vas  measured  are  In  the  sectorial  boxes  shorn  in 
Fig.  9.2. 

We  regard  the  conparlsons  of  theory  and  experiment  in 
Fig.  9.2  as  exploratory  and  not  definitive.  It  is  of 
course  striking  that  the  experiments  of  Metzner,  Uebler 
and  Fong  (1969)  do  appear  to  involve  a  vortlcity  of 
changing  type.  It  would  be  Interesting  to  6ee  if  thi6 
striking  type  of  experimental  result  could  be  repeated 
by  other  investigators  using  different  fluids  and 
experimental  arrangements.  We  hasten  to  add  that  the 
Separan  solution  used  in  the  experiment  is  not  an  Oldroyd 
model  and  surely  cannot  be  characterized  by  a  viscosity 
and  relaxation  time.  In  fact  only  special  models  give 
the  vortlcity  precisely  a6  the  quantity  which  changes 
type.  Ve  have  already  remarked  that  models  with  true 
viscosity;  e.g.,  retardation  times,  will  smooth  discon¬ 
tinuities,  with  only  a  little  smoothing  if  the  retarda¬ 
tion  "viscosity"  parameter  is  small.  Probably  all  the 
polymer  solutions  used  in  experiments  have  some  small 
smoothing.  In  viev  of  all  these  uncertainties  in  theory 
and  experiments  it  would  be  premature  to  make  strong 
claims . 
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10.  CHANNEL  FLOWS  WITH  WAVY  WALLS 

Yoo  and  Joseph  0984)  solved  the  problem  of  flow  of  an  upper  con- 
vected  Maxwell  model  through  a  channel  with  wavy  walls  linearized  for 
small  waviness.  It  is  easy  to  find  an  exact  solution  of  the  equations 
when  the  walls  are  flat.  The  velocity  profile  is  the  same  quadratic  one 
that  one  finds  in  Newtonian  fluids.  There  are  some  normal  stresses 
which  are  absent  in  the  Newtonian  case.  This  Poiseullle  flow  is  per¬ 
turbed  and  linearized.  The  characteristics  and  regions  of  hyperbolicity 
can  be  computed  without  specifying  the  nature  of  the  perturbation.  Yoo 
and  Joseph  went  further.  They  used  a  specific  perturbation,  the  ampli¬ 
tude  of  the  waviness  of  the  walls.  They  defined  the  linearized  problem 
for  small  waviness  and  solved  it  numerically.  The  Yoo-Joseph  paper 
gives  the  first  actual  computation  of  a  flow  with  change  of  type. 

We  are  going  to  outline  the  analysis  and  some  of  the  main  results 
of  Yoo  and  Joseph.  We  shall  express  the  equations  in  terms  of  the 
"Mach"  number  M  »  U/C  and  the  elasticity  number  E  -  nl/pd2  (this  is  in¬ 
dependent  of  U).  The  Wei3senberg  number  U0I/d  is  given  by  W  -  M/E. 


uu  ♦  vu  +  4i(p  "0  -  t  )  -  0, 
x  y  M2  Kx  x  y  ’ 

UV  ♦  VV  +  T7?(p  ■  t  “  Y  )  -  0, 
x  y  M2  Ky  x  y  ’ 

uo  *  vav  *  2( o*1  )u  -  2 tu  - - — , 

x  y  X  y  M/E 

UTv  *  V  Tw  "  +  ~  (0  +  1)v  - - — , 

X  Y  y  X  M/E 

uY  ♦  vY  -  2  tv  -  2 ( Y ♦  1 ) v  ■ — . 

X  Y  y  M/E 


(10.1a) 

(10.1b) 

(10.1c) 

(10. Id) 

(10. 1e) 

(10. If) 


We  shall  seek  and  find  a  solution  of  these  equations  satisfying  no  slip 
conditions  at  the  walls 
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u  ■  v  ■  0  at  y  ■  i(1  ♦  e  sin  nx)  (10.2) 

with  symmetric  streamlines 

|^-v-0aty-0.  (10.3) 

3y 

It  is  noteworthy  that  our  solution  is  completely  determined  by  data 
(10.2)  and  (10.3)  on  the  velocity  alone.  It  is  not  necessary,  and  it 
would  be  wrong  to  prescribe  more  about  velocity  or  stresses.  The  vorti- 
city  jo  -  curl  u  is  related  to  the  streamfunction 

JS  “  curl  H  "  «2*  ’  -[$  *  0]®2  (10-4) 

where  u(x,y)  satisfies 


(MV-o-I)u  ♦  2(M*uv-t)(ii  ♦  (M*v*-Y“1  )<*) 

xx  xy  yy 


♦ 

-p  ♦  JLU 

0)  ♦ 

-p  ♦  JLy 

X  /  E 

X 

y  /E 

♦  (u  -v  )(o  +  Y  ♦!  +  T  )  ♦  (u  *v  )( Y  -o  ) 

x  y  xy  xy  xx  yy  y  x  yy  xx 


♦  o  u  ♦  2t  u  ♦Yu  -ov  -2xv  -Yv  -0. 
y  xx  y  xy  y  yy  x  xx  x  xy  x  yy 


(10.5) 


Equation  (10.5)  may  be  written  as 

Au  ♦  2Bui  ♦  Cu»  ♦  L.O.T.  -  0  (10.6) 

xx  xy  yy 

where  the  terms  L.O.T.  are  of  lower  order  for  hyperbolic  analysis  (see 
JRS,  1 98 i* ) .  Characteristic  directions 
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1 1 
dx 


A  -  M*u*  -  o  -  1 , 


B  -  M^JV  -  t , 


C  -  M*v*  -7-1 


for  the  vorticity  exist  whenever 


(10.7) 


def 

I  «  t2  -  2M2tuv  -  (1  +  Y)(1  +  o)  ♦  M*v*(1  ♦o)  ♦  MV(H-Y)  >  0.  (10.8) 

The  expression  (10.8)  is  expressed  in  terms  of  unknown  velocity  and 
stress  fields.  The  criterion  (10.8)  for  hyperbolicity  can  be  satisfied  in 
some  regions  of  flow  and  not  in  others.  The  border  £  »•  0  between  the 
elliptic  and  hyperbolic  regions  of  flow  is  like  the  sonic  line  in  gas 
dynamics.  Across  this  line  the  equations  are  said  to  change  type. 

Equation  (10.1)  and  all  the  equations  of  this  section  are  general  in 
that  they  apply  to  every  plane  problem,  not  just  the  channel  flow  prob¬ 
lem  Introduced  in  Section  1. 

Now  we  shall  solve  the  governing  equations  for  flow  in  a  channel 
with  straight  walls  e  -  0. 


(u0,v0)  -  (l-y2,0), 


(Po.t«)  -  -2M(x,y)/E, 


vo0,y0)  -  (8M2Ey2,0), 

-  2y  (10.9) 


T.ne  basic  motion  depends  exclusively  on  the  Weissenberg  number  W  »  M/E 
measuring  the  size  of  stresses.  The  solution  (10.9)  is  relatively  fea¬ 
tureless  and,  in  particular,  it  gives  no  indication  of  hyperbolicity. 

Now  we  consider  any  plane  perturbation  of  (10.9).  The  problem  with 
wavy  walls  is  one  such  perturbation,  but  there  are  infinitely  many 
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others.  We  may  linearize  the  formula  (10.3)  for  the  characteristics  of 
any  flow  slightly  perturbing  the  Poiseuille  flow  (10.9).  The  charac¬ 
teristics  for  all  these  perturbations 


dy  1 
dx  “  M 


* 


-2y/£  ±  ■  4Ey2  "  "m? 


-(y2-1)2  ♦  8EyJ  +  jp 


(10.10) 


are  defined  in  terms  of  quantities  defined  for  the  basic  flow  (10.9)  and 
are  given  once  and  for  all,  independent  of  the  perturbation. 

Equation  (10.10)  shows  that  flows  perturbing  plane  Poiseuille  flow 
can  exhibit  a  change  of  type  with  a  "sonic"  line  £  -  0  given  by 


Ky2)  -  (y2-i)2  -  2Ey2  -  jj,  -  o  (-1  s  y  s  i). 


Since  E(y2)  is  monotonically  decreasing,  it  has  a  maximum  at  y2  *  0  and 


E(°)  -  1  -  jp  >  0 


if  and  only  if  the  viscoelastic  "Mach"  number  M  >  1.  The  "sonic"  line 
across  which  the  flow  changes  type  is  y  •  y*  where  E(y*2)  -  0 


1  ♦  2E 


-  2^EJ 


The  linearized  problem  for  small  e  is 

u  +  v  -  0, 
x  y 

(1-y2)ux  -  2yv  +  fp(Px~ °x" V  ”  °* 

( 1  ~y2) V  ♦  77j(p  -t  -v  )  -  0, 

J  x  M2  *y  x  y 

(1-y2)ox  ♦  l6M2Eyv  -  2(8M  JEy J+ 1  )u^  ♦  ^  y  t  ♦  ^Myu^/ E* o/M/ E 


(10.11a) 

(10.11b) 


(10.11c) 


0, 

(10. lid) 
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d-yJ)Yx  ♦  4Myvx/E  -  2vy  ♦  Y/M/E  -  0,  (lO.lle) 

(1-y*)x  -  2Mv/E  ♦  2y Y  -  (8M2Ey2*1)v  -  u  ♦  t/M/E  -  0  (10. Ilf) 

x  x  y 


where 


I 


u(x,±1)  -  2  3in  nx 
v(x,±1)  -  0, 


v(xf0)  -  uy(x,0)  -  0.  (10.12) 

Yoo  and  Joseph  solved  the  equations  (10.11)  using  only  the  velocity  data 
(10.12)  and  the  method  of  separation  to  reduce  these  equations  to  ordi¬ 
nary  differential  equations  in  y.  The  ordinary  differential  equations 
were  solved  numerically.  The  vorticity 


def 

u  "  vx  '  Uy  *  ~V2’1' 
satisfies  an  equation 


(10.13) 


[(.-/■)■  -  9Ey*  -  £].xx  •  4^  y»xy 


— ~(1  -y2) 
M/E 


v 

M/E 


l6Eyu  ♦ 
xx 


4/E 
M  Uxy 


2(1-y2)vx 


I 


I 


» 


♦ 


2yY 


yy 


2yo  ♦  4Y  ] 
xx  y 


(10.14) 


of  changing  type  whose  characteristics  are  given  by  (10.7). 

In  Figs.  10.1,  10.2  and  10.3  we  have  graphed  vorticity  and  stream¬ 
lines  for  M  -  10.  In  10.1a  and  10.2a  we  superimposed  characteristics 
(light)  on  zero  vorticity  lines.  We  first  observe  that  when  W  is  small 
the  decay  of  the  vorticity  is  rapid.  When  W  is  large  the  vorticity 
decays  only  very  slowly  with  oscillations  all  the  way  to  the  center. 
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2.  FLOW  OVER  A  FLAT  PLATE 

The  problem  of  flow  around  a  flat  plate  could  be  thought  to  be 
loser  than  flow  around  a  body  to  a  perturbation  of  uniform  flow.  In 
act,  the  theory  of  slender  bodies  in  aerodynamics  is  a  perturbation  of 
niform  flow  which  is  perturbed  less  by  slender  bodies  than  by  fat  ones, 
he  aerodynamic  theory  works  well,  but  only  because  of  the  flow  follow- 
ng  condition  in  which  the  fluid  is  required  to  slip  along  the  slender 
ody.  In  this  chapter  I  am  going  to  give  an  exact  supercritical  theory, 
alid  for  all  fluids  with  instantaneous  elasticity,  such  that  the  vorti- 
ity  is  prescribed  arbitrarily  along  the  flat  plate,  with  a  vanishing 
orraal  component  of  velocity  on  the  along  the  flat  plate,  with  a  vanishi¬ 
ng  normal  component  of  velocity  on  the  plate  surface.  Presumably  the 
o  slip  condition  on  the  plate  can  also  be  satisfied  by  choosing  the  cor- 
ect  vorticlty  distribution  on  the  plate.  In  the  nonlinear  problem  there 
ill  be  a  diffusive  subcritical  region  near  the  plate  and  a  "transonic" 
urface  outside  where  the  viscoelastic  Mach  number  passes  through  one. 
his  type  of  important  problem  has  yet  to  be  considered  and  solved  in 
he  theory  of  flow  of  viscoelastic  fluids. 

We  now  consider  two  dimensional  flow  past  a  flat  plate.  The  plate 
3  on  the  half  line  y  -  0,  x  Z  0.  The  velocity  components  corresponding 
o  (x,y)  are  (u,v)  and  there  is  only  one  component 


3v  3u 
**  3x  9y 

f  vorticity.  The  velocity  components 


(12.1) 


(12.2) 


ay  be  obtained  from  a  stream  function  Y,.  We  suppose  that  the  velocity 
i,v)  vanishes  at  infinity  and  that 


(u,v)  «  (-11,0)  on  y  =*  0,  x  2  0. 


(12.3) 


The  vorticity  of  the  perturbed  flow  satisfies 


w 
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ZERO 

VORTICITY 

M  U/C, 

C  -  /G(0)/p 


Fig.  11.1  Leading  characteristics  for  vorticity  in  a  plane 
uniform  flow. 


TERMINAL 

VELOCITY 


Fig.  11.2  The  jump  in  the  terminal  velocity  for  an  air 
bubble  rising  in  a  liquid.  The  jump  can  be 
large,  say  a  five-fold  increase. 


where  d  la  a  constant  and  U  Is  In  the  direction  x.  The  characteristics 
form  a  net  of  straight  lines. 

If  we  suppose  that  a  small  two  dimensional  body  perturbs  uniform 
flow  when  U  >  C  then  there  would  be  an  undisturbed  region  in  front  of 
the  body  which  could  not  be  reached  by  disturbances  of  the  vorticity 
traveling  at  velocity  C.  In  a  linear  theory  the  first  changes  in  the 
vorticity  would  occur  across  the  leading  characteristics  which  form  an 
angle  like  the  Mach  angle  of  gas  dynamics  (see  Fig.  11.1).  We  could  call 
the  undisturbed  region  of  uniform  flow  a  "region  of  silence". 

No  one  has  yet  solved  a  problem  of  supercritical  flow  over  a  body, 
even  in  the  linearized  case  and  no  one  has  looked  for  a  region  of  silence 
in  experiments. 

An  air  bubble  rising  in  a  liquid  will  reach  a  steady  terminal  veloc¬ 
ity.  The  larger  the  bubble  the  larger  is  the  terminal  speed.  Astarita 
and  Apuzzo  (1965)  were  the  first  to  notice  that  the  terminal  velocity  of 
the  bubble  Is  not  a  smooth  function  of  the  bubble  volume.  There  is  a 
critical  volume  which  is  associated  with  a  jump  in  the  terminal  velocity, 
as  in  Fig.  11.2. 

This  phenomenon  has  also  been  studied  by  Calderbank,  Johnson  A 
Loudon  (1970).  The  fluids  used  in  these  experiments  are  water-based 
polymers.  The  lower  value  of  the  critical  terminal  velocity  ranges  from 
0.1  to  10  cm  per  sec.  These  values  are  not  unlike  values  which  we  think 
are  typical  for  wave  speeds  in  these  polymers  (see  Section  4).  The  rea¬ 
sons  for  the  abrupt  rise  in  the  terminal  velocity  are  not  understood 
Maybe  the  abrupt  change  is  associated  with  a  change  of  type.  Then  we 
might  expect  to  see  an  abrupt  change  near  M  -  1.  For  Maxwell  models  M 
-  RW  where  R  is  the  Reynolds  number  and  W  is  the  Weissenberg  number. 
/Ian a  and  Leal  (1978)  give  values  of  RW  which  vary  by  decades  at  critical¬ 
ity.  Their  method  of  calculating  was  not  clear.  They  seemed  to  have 
used  some  normal  stress  data  which  could  give  very  inaccurate  values  for 
W.  if  they  are  nearly  right  it  will  be  hard  to  support  the  idea  that  the 
ibrupt  rise  is  associated  with  a  change  of  type.  The  problem  is  open. 
To  settle  this  problem  we  need  to  have  accurate  values  for  the  wave 
velocity  used  in  experiments. 
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pU* 


3x 


1 

U7 


£(l)  dx 


♦ 


’2* 


£<*)  dx 


(11.12) 


The  analysis  of  supercritical  flow  is  conveniently  framed  for  equa¬ 
tions  of  second  order.  We  differentiate  (2.12)  with  respect  to  x  and 
find  that 


(M*-1 ) 


G"(0)  3I 

to*  ~  uctoT  3x 


Vi  ♦ 


G"(0)  r 
U*G(0)  * 


•  uck  U0!2?]  V  *  V  a»- 


(11.13) 


This  equation  (11.1 3)  changes  type  from  elliptic  to  hyperbolic  when 
M*  -  pUVG(O)  -  U*/C2 


Increases  through  one. 

For  Maxwell  models  (11.13)  reduces,  using  (11. *0,  to 


this  equation  changes  type  when  M*  -  U*/C*,  C*  -  n/pX  increases  through 
one. 

The  vorticlty  of  steady  flows  with  instantaneous  elasticity  perturb¬ 
ing  uniform  flow  will  change  type  from  elliptic  to  hyperbolic  when  the 
ratio  of  the  velocity  of  the  free  stream  to  the  velocity  of  propagation 
of  shear  waves 

into  a  fin’d  at  rest  exceeds  unity.  Characteristics  for  vorticlty  in 
plane  flow  are  given  by 


'v 


_• 


M 
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♦  u 


G(t-T)A£(jt,  T)dx. 
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(11.7) 


The  x  derivatives  under  the  Integral  may  be  replaced  with  x  derivatives. 
Equation  (11.7)  shows  that  states  of  uniform  vorticity  are  solutions  of 
the  equations  of  motion  of  fluids  with  instantaneous  elasticity  which 
perturb  uniform  flow.  Potential  flow  is  possible. 

For  steady  flows  it  is  convenient  to  change  variables  in  the  inte¬ 
grals.  We  write 

X  -  x  -  U(t--r),  dx  *  Udx  (11.8) 

Equation ,(1 1 .3)  may  be  written  as 


G(0)Aj[u(x)] 


Atu(j^)]  dx 


(11.9) 


Equation  (11.6)  is  now 


3i 

pU  —  -  curl  div  x 


(11.10) 


and  (11.7)  becomes 


pU* 


(ii.ii) 


which,  after  integrating  the  first  term  under  the  integral  by  parts, 
becomes 
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[t 

3t  3t 

t=  ♦  U  t=  -  G(0)A[u(x,t)]  ♦  G'(t-T)A[u(j^,  t)]  d t 

at  aX  —  —  — 


(11.3) 


where  the  last  term  is  of  lower  differential  order  for  hyperbolic  analy¬ 
sis. 

For  all  the  Oldroyd  models  the  perturbation  of  uniform  motion  leads 
to  the  linearized  Maxell  model  with 


G(s)  -  7  exp(-s/X) 


(11.4) 


The  rate  equation  for  a  Maxwell  model  in  a  flow  perturbing  uniform  flow 


3t  3t 


f  *"5'?  ACu(x,t)]  -  =. 


(11.5) 


The  equations  of  motion  are 


f  3u  3ul 

(a!  *  u  taj  '  '7p  *  dIV  1- 


The  vorticity  £  -  curl  u  satisfies 


{I  * u  I] 


curl  div  t. 


(11.6) 


Using  (2.1)  and  curl  div  A[u(x,  t)]  -  A£(^,  t) 


where 


V.  7,2 


JL  ♦  * 

3p  3P’ 


we  get 


V 
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ty  to  the  wave  speed  is  greater  than  one,  and  the  wave  velocity  which  is 
a  property  of  the  fluid,  cannot  be  controlled.  It  follows  that  even  if 
the  flow  away  from  the  body  is  supercritical  the  flow  near  the  body 
will  always  be  subcritical,  and  the  underlying  problem  is  really  one  of 
"transonic"  type. 

We  are  going  to  show  that  the  equations  which  govern  perturbations 
of  uniform  flow  have  a  simple  form,  depending  on  the  density  p  and 
relaxation  modulus  G(s)  alone,  and  are  model  independent.  The  flow  near 
the  body  depends  on  the  constitutive  model. 

The  stress  in  a  simple  incompressible  fluid  may  be  decomposed  into 
an  isotropic  and  a  determinate  part.  The  determinate  part  is  the  consti¬ 
tutive  equation  which  relates  stress  and  deformation.  In  linearization 
at  uniform  motion,  the  determinate  stress  in  fluids  with  instantaneous 
elasticity  is  given  by  an  integral  (see  JRS,  1984) 


T 


t 

G(t-t)A[u(j£,  t)]  dt 


(11.1) 


where 


x  -  U(t-x) 

y 

z  J’ 

,  T  3ui  3ui 

Atu]  .  Vu  .  VHT,  jy  .  _  . 

and  G  is  smooth,  positive  and  monotone  decreasing.  The  speed  of  various 
kinds  of  shear  waves  (Coleman  and  Gurtin,  1968,  Narain  and  Joseph,  1982) 
into  a  fluid  of  density  p  at  rest  is 


C 


(11.2) 


The  "Mach"  number  is  U/C.  By  differentiating  (2.1)  with  respect  to  t, 
holding  x  fixed,  we  find  that 


A  i 

M  -  ^  -  U// n/ A  p 

exceeds  one  (see  Eq.  ( 2 - ^ .  Here  U  is  the  velocity  of  the  unperturbed 
uniform  flow,  and  C  Is  the  wave  velocity  for  propagation  of  shear  waves 

In  a  Maxwell  fluid.  Ultman  and  Denn  (1970)  did  not  notice  that  It  Is 

precisely  the  vorticity  which  changes  type.  They  refer  to  their  lineari¬ 
zation  as  an  Oseen  approximation.  Oseen  introduced  his  linearization 

around  uniform  flow  for  slow  viscous  flow,  because  Stokes  equations  have 
no  solution  for  flow  around  bodies  in  two  dimensions.  Oseen's  equations 
do  not  change  type.  Ultman  and  Denn  (1971)  use  Oseens  approximate 
method  to  compute  subcritical  flow  and  they  say  that  their  calculations 
agree  with  their  experiments  when  the  fluid  parameters  are  properly 
chosen. 

In  a  second  paper,  Ultman  and  Denn  (1970)  consider  the  supercritical 
flow  but  they  do  not  give  experimental  results  and  they  do  not  discuss 
or  try  to  solve  their  equations.  They  attempt  to  correlate  some  experi¬ 
mental  observations  of  D.  F.  James  (1967)  with  the  change  of  type. 
James  observes  a  sudden  change  in  the  slope  of  the  heat  transfer  curve 
as  a  function  of  velocity.  This  happens  at  a  critical  velocity  which  Tor 
the  Polyox  solution  used  by  James,  was  about  1  cm/sec.  It  is  not  clear 
from  the  graphs  how  abrupt  this  change  of  slope  is,  but  there  is  a 
change  of  slope.  Ultman  and  Denn  (1970)  also  suggest  that  the  transi¬ 
tion  from  subcritical  to  supercritical  flow  might  explain  abrupt  changes 
in  the  drag  coefficient  they  say  was  observed  by  A.  Fabula  (1966). 
Again,  the  idea  is  that  the  critical  velocity  at  transition  is  the  wave 
speed  C.  They  make  an  estimate  of  C  from  a  molecular  theory  and  corre¬ 
late  this  prediction  with  the  data  of  James.  Of  course,  any  such  esti¬ 
mate  can  at  best  be  expected  to  give  an  order  of  magnitude,  since  the 
molecular  theory  is  coarse  and  the  fluids  used  in  the  experiment  are  not 
Maxwell  fluids. 

Recently  Ambari,  Deslouis  and  Tlbollet  (198*4)  have  considered  mass 
transfer  in  the  flow  of  a  viscoelastic  fluid  around  a  cylinder.  They 
also  find  a  critical  value  at  which  the  mass  transfer  undergoes  an  ano- 
molous  transition. 

Oseen's  methods  work  because  3low  flow  around  bodies  is  in  some 
sense  a  perturbation  of  uniform  motion.  The  case  of  fast  flow  is  dif¬ 
ferent  because  the  uniform  stream  is  not  small,  but  the  velocity  on  the 
body  must  vanish.  In  supercritical  flow  the  ratio  of  the  stream  veloci- 
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Yoo  and  Joseph  did  an  asymptotic  analysis  which  give  the  frequency  of 
the  oscillations  and  other  properties  of  the  solution.  The  second  prop¬ 
erty  is  that  regions  of  positive  and  negative  vorticity  which  are  not 
damped  are  swept  out  along  characteristics.  It  would  be  hard  to  under¬ 
stand  these  solutions  from  the  streamline  plots. 

The  same  type  of  hyperbolic  dynamics  which  Yoo  and  Joseph  found  in 
channels  will  occur  in  pipes.  Mark  Ahrens  is  working  on  this  problem. 
The  change  of  type  which  occurs  in  the  center  of  a  pipe  may  have  appli¬ 
cations  in  the  problem  of  delayed  die  swell.  The  phenomenon  of  delayed 
die  swell  is  not  well  known.  At  low  speeds  the  Jet  will  spread  near  the 
exit  of  the  Jet,  as  in  Fig.  10.H(a)  and  10.5(a).  At  yet  higher  speeds  the 
Jet  doe3  not  spread  near  the  exit,  the  swell  is  delayed,  as  in  Fig. 

10. U(b-d).  Fig.  10.5  shows  a  form  of  delayed  die  swell  which  we  have 
seen  repeatedly  in  our  own  experiments.  The  delayed  swell  seems  to 
occur  at  a  critical  speed,  not  so  different  than  what  one  might  expect 
from  a  change  of  type.  Of  course,  the  reason  for  the  delayed  swell  is 
not  understood.  The  form  of  the  jet  reminds  one  of  a  hydraulic  jump, 
which  is  the  shock  phenomenon  corresponding  to  shocks  in  gas  dynamics. 

11.  PROBLEMS  ASSOCIATED  WITH  THE  FLOW  OF  VISCOELASTIC  FLUIDS 
AROUND  BODIES 

It  is  perhaps  not  unreasonable  to  think  that  far  from  the  body  we 
have  only  a  small  perturbation  of  uniform  motion  with  constant  velocity 
Ue^  in  the  direction  x.  The  body  is  stationary. 

The  vorticity  of  all  steady  flows  of  viscoelastic  fluids  with  in¬ 
stantaneous  elasticity  which  perturb  uniform  flow  can  change  type. 

To  be  precise,  the  linearized  equations  for  the  vorticity  of  flows 
perturbing  uniform  flow  vanish  on  states  of  constant  vorticity  and  the 
type  of  this  equation  in  steady  flow  changes  when  the  ratio  of  the  free 
stream  velocity  of  the  wave  3peed  increases  through  unity.  This  criter¬ 
ion  for  a  change  in  type  may  be  expressed  in  terms  of  a  viscoelastic 
Mach  number. 

Ultman  and  Denn  (1970,  1971)  consider  the  equations  for  two-dimen¬ 
sional  steady  flow  of  an  upper  convected  Maxwell  fluid.  They  linearize 
at  a  motion  with  uniform  velocity  and  zero  stress,  and  they  show  that 
these  linearized  equations  change  type  when  a  viscoelastic  "Mach"  number 
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Figure  10.4  (After  Brenschede  and  Klein,  1970).  Delayed  die  swell 
in  polvisobutilene  solution  is  toluol. 

D  -  4Q/ irR1  where  5  is  the  volume  flow  rate  and  R  is 
the  pipe  radius.  The  critical  value  D  is  14  x  10- 
reciprocal  seconds. 

(a)  D  =  8  x  10-  (b)  D  =  =  15.2  x  10- 

(c)  D  -  22.8  x  10-  (d)  D  -  29  x  1 0H 


i 


-4 


®  Figure  10.5  (After  H.  Giesekus,  1968).  Delayed  die  swell  in  a 

.5i  polyacrylamide  solution  in  water. 

(a)  slow  speed,  (b)  post  critical  speed, 

(o )  high  sp*’ed 


For  the  Maxwell  model  (11.4)  we  get  a  telegraph  equation 


We  may  write  this  telegraph  equation  in  dimensionless  form 
ili  +  M2  l£  1  j£c  n 

3?  +  B*W  to  '  ?*  3^  “  0 
where  (x,y)  -  (x/lty/£),  l  is  a  scale  length 

8*  -  Ml  -  1  >  0 

and 

W  -  UX/fc 


(12.5) 


is  a  Weissenberg  number  based  on  l. 

Equation  (12.5)  can  be  approximated  by  a  wave  equation  without 
damping  when  W  ■*  «  and  M  >  1  is  fixed.  For  this  same  M  we  may  expect 
rapid  damping  of  the  backward  heat  equation  when  W  is  very  small. 

The  other  prescribed  conditions  are  that  the  vorticity  must  vanish 
far  from  the  plate  and  since  v  *  0  on  the  plate, 

C  *  “J~  on  y  -  0,  x  2  0  (12.6) 

We  need  vorticity  fields  satisfying  (12.4)  or  (12.5)  over  R*  and  (12.6)  and 
(12.3)  on  the  plate. 

The  stream  function  and  vorticity  are  related  by  a  second  order 
equation 


v  V,"  ^  .•  ,■ 11  ■  ■  ■  9>m?  "  ^  V  "T  >  ■>  • ..‘V--  ”■ J-. 


-c(y.x). 
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(12.7) 


We  may  therefore  expect  that  corresponding  to  any  good  field  c(x,y)  we 
may  solve  (12.7)  subject  to  the  conditions  that  f  vanishes  far  from  the 
plate  and 

v  ■  I-  •  0  on  y  ■  0,  x  J  0.  (12.8) 


In  general,  we  cannot  expect  that  for  any  field  t(y,x) 

II  ■  ^  on  y  ■  0,  x  i  0.  (12.9) 

It  follows  that  the  existence  of  a  solution  of  the  linearized  flat  plate 
problem,  in  both  the  subcritical  M  <  1  and  supercritical  M  >  1  case  de¬ 
pends  on  finding  a  special  prescription  of  the  vorticity  c  on  the  plate 
(12.6)  such  that  (12.9)  is  verified. 

I  wish  now  to  consider  the  supercritical  case  and  to  indicate  the 
method  by  which  our  linearized  problem  may  be  resolved.  It  will  be  con¬ 
venient  to  start  this  discussion  with  an  analysis  of  the  telegraph  equa¬ 
tion  (12.5)  for  M  >  1.  To  put  this  problem  into  canonical  form  we 
change  variables  setting 


and  find  that  c(z,t)  -  c(y,x)  satisfies 


(12.10) 


£L  ♦ 

at1  at  P 


o 


(12.11) 


We  may  expect  that  this  hyperbolic  problem  will  not  allow  the  plate  to 
influence  the  flow  upstream.  We  therefore  seek  solutions  £(z,t)  such 
that  the  upstream  vorticity  vanishes; 

t(z,t)  -  0  for  t  <  0,  -•  <  z  <  - 


This  implies  that  "initially" 


C(Z|0)  -  "  0 


(12.12) 


whereas 

f  • 

I 

£(0,t)  -  g(t)  (12.13) 

is  prescribed  for  t  >  0.  Asymptotically,  for  large  |z|  we  require  that 

I  l  *  0. 

The  problem  (12.11),  (12.12)  and  (12.13)  can  be  solved  by  Laplace 
transform  techniques.  A  well  known  solution  of  this  problem  with  uni¬ 
form  vorticity  on  the  plate 


f(t)  -  H(t)  - 


0,t<0 

1,t>0 


(12.00 


has  been  given  by  Carlslaw  and  Jaeger  (1963)  in  the  form 


t(z.t) 


I 


I,  (^/  o4-za)  do 


H(t-z) 


(12.15) 


where  I,  is  the  modified  Bessel  function  of  the  first  kind.  The  vorti¬ 
city  £(-z,t)  on  the  bottom  of  the  plate  is  given  by  02.15)  with  -z  re¬ 
placing  z.  This  solution  for  flow  over  a  flat  plate  was  given  by  B. 
Caswell  (1976)  in  a  3tudy  of  the  effects  of  a  leading  edge  singularity. 

We  next  note  that  the  field  c(y,x)  -  c(z,t)  can  now  be  inserted  into 
(12.7),  which  is  an  elliptic  problem  leading  to  a  nearly  everywhere  dif¬ 
ferentiable  T,  even  when  C  has  simple  discontinuities,  as  in  Figure  1 . 
We  solve  (12.7)  subject  to  02.8),  but  the  solution  will  not  satisfy 
(12.9).  We  can  hope  to  satisfy  (12.8)  and  (12.9)  simultaneously  by  pre¬ 
scribing  the  perfect  vorticity  distribution  f(t)  on  the  plate. 

To  get  the  £(t,z,g)  corresponding  to  different  prescriptions  (12.13) 
of  the  vorticity  f(t)  on  the  plate,  we  could  use  the  method  of  Duhamel 
type  integrals  introduced  for  start  up  problems  by  Narain  and  Joseph 
0983).  For  these  integrals  we  need  to  superpose  using  the  fundamental 
singular  solution  of  (12.11),  (12.12)  and  (12.13)  when 


g(t)  -  6(t-t) 


(12.16) 


is  a  Dirac  function.  It  is  easy  to  see  and  not  hard  to  prove  that  the 
required  solution  is  the  time  derivative  of  the  step  function  problem 
just  derived  when 


C(0,t)  - 


1  for  t  >  t 
0  for  t  i  t' 


Using  the  aforementioned  method  I  find  that  the  solution  of 
(12,12)  and  (12.13)  Is 


C(t,z) 


t-z 

g(t)  |£(z,t-T)dx  ♦  e  z/2  g(t-z) 


where 


f(z,t)  -  e 


-z/2 


-o/2 


I|(  jj/  o*-z5) 
/ a2-zi 


do 


The  solution  c(y.x)of 

14  ♦  Ml  li  34 

3x5  rw  3x  B1  By7 


0  . 


q(0,x)  -  g(x), 
c(y.o)  -  0 


is 


rx-By 


<;(  y.x) 


g(r)  |^(y,x-T)dx  ♦ 


ISs] 


g(x-By) 


(12.17) 

(12.11), 

(12.18) 


(12.19) 

02.20) 


where 
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The  amplitude  of  ;  on  *  ■  ly  la  given  by  exPt“2^J^*  The  aBPlltu<le 
decays  rapidly  when  W  is  small  (see  Fig.  12.1). 

Turning  next  to  the  stream  function,  we  find  that  the  solution  of 
(12.7)  with  c(y,x)  given  by  (12.20)  and 

i|»(x,0)  -  0  for  x  >  0  (12.21) 

is 

♦  ■  “ j  |  G(y,y,|x,x0)c(y0,x0)dy0  dx,  (12.22) 


where  R_  is  the  complement  to  the  line  x  >  0,  y  -  0.  The  Green  function 
for  this  domain,  vanishing  on  the  plate  is  given  by  (10.1.40)  on  p.  1208 
of  Morse  &  Feshbach  (1953).  Me  have  to  divide  their  solution  by  4w. 

G(y,y0|x,x0)  -  G(r,r0|  $,$,) 


i  i  i  •-[=?]  -MC- 


r  <  r0 


1 1  i  4?] 


r  >  r« 


(12.23) 


where  (x,y)  -  r(cos$,sin$) 

The  solution  (12.7),  satisfying  (12.21)  is  given  by  (12.22)  with 
C(y.x)  given  by  (2.20).  The  fluid  will  not  slip  on  the  plate  if 

U  -*^(x,0)  when  x  >  0;  that  is, 

oy 
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-u  -  ^  y-o  G^*y«lx»x«^<y®*x«)dy« dx®  02.2*0 

/  Jr.  |x>0 

We  have  to  choose  the  prescribed  plate  vortlcity  g(x)  to  satisfy  (12.24). 
I  have  no  guarantee  that  there  is  a  g(x)  satisfying  (12.24).  Edmond 
O'Donovan  is  trying  to  find  a  numerical  approximation  for  the  g(x)  which 
satisfies  (12.24). 

It  is  probable  that  the  approach  which  I  have  taken  in  this  problem 
can  be  generalized  to  fluids  satisfying  (12.4).  This  equation  applies  to 
all  fluids  with  instantaneous  elasticity.  Such  a  study  would  show  how 
the  flow  depends  on  the  kernel  G(s).  The  Laplace  transform  methods  used 
by  Narain  and  Joseph  (1982,  1983)  to  study  the  linearized  dynamics  of 
shearing  motions  perturbing  rest  in  viscoelastic  fluids  are  appropriate 
for  this  study  too.  As  in  their  work  I  get  the  following  problem  for 
the  transform 


c(y.u)  -  I  eUX  c(x.y)  dx 
Jo 

where  c(y.«)  vanishes  for  large  y  and  satisfies 

k(u));  -  g(u>) 


(12.25) 


(12.26) 


and 


kU)  -  «XpU*  -  G(0) J  - 


(0  °"W ds- 


e"wS  G 


M  *• 


One  finds  that 


Ku,  .  SM  ♦  o 

U  I 


.  G(O)  .  G*(0)  .  nf  1  1 

8<«>  ■  —  *  i>r  °L?J 

if  G"  is  integrable.  The  character  of  the  solutions  is  determined  by  the 
symbol  of  the  operator.  We  should  therefore  look  at  the  problem  which 
arises  at  large  u 


h-  -  «•»  -  °-¥\i  -  [*?  *  w]  $ 


[8V  ‘  o^T  •  0 


which  is  the  transform  of  the  equation 


■  a,  .  P»9'(0?  ,  .  , 

8  Cxx  G*(0)  Cx  Cxx’ 


(12.27) 


From  this  equation  we  expect  that  the  solution  of  the  present  problem  is 
a  simple  wave  as  in  Figure  12.1  with  c  ■  0  to  the  left  of  the  line 
x  -  By  on  the  top  of  the  plate,  and  x  -  -By  on  the  bottom  plate.  The 
amplitude  of  the  wave  on  x  •  By  is  given  by 


,„_fxpUG'(0)1 

6XPL  82G*(0)  J 


60 


This  amplitude  tends  to  zero  as  G'(0)  -*  When  G(0)  >  0  and  G'(0)  -  -• 

we  have  the  case  of  singular  kernels.  Renardy  (1982)  showed  that  for 
some  special  singular  kernels  shear  waves  would  propagate  into  a  fluid 

at  rest  with  the  wave  velocity  C  -  /G(0)/p,  with  a  zero  amplitude  at  the 
front  and  C*  connection  with  non-zero  solution  across  the  front.  These 
properties,  as  well  as  the  analytic  smoothing  of  sharp  fronts  with  small 
viscous  terms  (small  retardation  times)  occur  in  the  present  steady 
problem  as  well  as  in  the  theory  of  propagation  of  waves. 

Supercritical  flow  past  a  flat  plate  cannot  be  considered  to  be  a 
perturbation  of  uniform  flow.  If  it  is  supercritical,  the  free  stream 
velocity  is  finite,  perhaps  large,  but  the  fluid  must  come  to  rest  on  the 
plate.  Near  the  plate,  the  velocity  will  be  small  and  the  local  "Mach" 

number  less  than  one  even  when  M*  -  pUa/G(0)  >  1.  The  governing  prob¬ 
lem  near  the  plate  is  therefore  elliptic,  or  at  least  not  hyperbolic.  It 

may  be  true  that  when  M  is  large  the  solution  of  the  nonlinear  problem 
is  close  to  the  linear  solution  (if  it  exists)  except  in  regions  -immedi¬ 
ately  near  the  plate.  This  might  suggest  boundary  layers,  but  that 
thought  should  be  eschewed.  We  should  instead  think  of  a  narrow  sub- 
critical  region  near  the  plate  which  goes  supercritical  at  the  "sonic" 
line  which  could  also  be  close  to  the  plate  when  M  >  >  1.  We  have 

therefore  to  consider  a  "transonic"  type  of  problem  and  not  a  boundary 

layer  type  of  problem.  We  know  almost  nothing  about  such  problems  for 
viscoelastic  fluids. 

13.  NONLINEAR  WAVE  PROPAGATION  AND  SHOCKS 

It  is  well  to  motivate  this  chapter  by  reminding  the  reader  of  the 
huge  differences  between  linear  and  nonlinear  theories  of  gas  dynamics. 
In  the  linearized  theory  there  are  no  essential  differences  between  rar¬ 
efaction  and  compression  waves.  These  waves  propagate  according  to  the 
wave  equation  without  change  of  form.  In  the  nonlinear  theory  an  impul¬ 
sive  rarefaction  will  be  smoothed  by  nonlinear  effects  and  a  smooth  com¬ 
pression  will  shock  up. 

Compression  waves  are  impossible  in  incompressible  fluids.  Instead 
we  may  perhaps  speak  of  waves  of  shear  or  of  waves  of  vorticlty.  Very 
little  is  known  about  the  nonlinear  effects  in  the  flow  of  viscoelastic 
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fluids.  We  should  like  to  know  the  answers  to  the  following  questions. 
We  are  given  a  constitutive  equation  in  the  class  with  instantaneous 
elasticity:  Different  subclasses  in  this  class  can  lead  to  different 

results.  We  want  some  classification  of  results. 

(i)  Suppose  that  we  are  given  smooth  data.  Is  there  a  shock  up? 

(li)  When  there  is  a  shock,  what  variables  (vortlclty,  velocity,  dis¬ 
placement)  become  discontinuous? 

(iii)  We  are  given  impulsive  data.  Is  it  possible  that  the  nonlinear 
terms  smooth  discontinuities? 

We  first  consider  the  problem  of  wave  propagation  in  fluids  under¬ 
going  rectilinear  shear  flow.  The  formulation,  due  to  Coleman  and  Noll, 
is  embodied  in  the  representations  for  shear  flow  shown  in  (1.2)  of  the 
Appendix,  by  M*  Slemrod  to  this  paper.  Coleman  and  Gurtin  (1968),  fol¬ 
lowing  earlier  work  on  longitudinal  acceleration  waves  in  compressible 
materials  (1965),  showed  that  the  amplitude  of  a  jump  discontinuity  in 
the  fluid  acceleration  satisfies  a  simple  nonlinear  differential  equation 
of  Bernoulli  type  with  coefficients  determined  by  the  instantaneous 
value  of  the  relaxation  kernel  at  the  wave  and  second  order  instantane¬ 
ous  modulus  evaluated  at  the  wave.  The  instantaneous  value  of  the 
relaxation  kernel  is  designated  as 

G(k.O) 

where  k  is  the  3hear  rate  at  the  wave  and  G(k,s)  arises  as  the  kernel  of 
the  integral  representation  implied  by  the  Riesz  theorem  for  first  func¬ 
tional  derivatives  of  the  functional  t  (•)  of  (1.2)  in  Appendix  A  evalu- 

a 

ated  in  the  weighted  L^CO,*]  spaces  of  Coleman  and  Noll.  This  kernel 
reduces  to  G(s)  when  the  wave  advances  into  a  region  at  rest.  The  speed 
of  the  wave 

c . 


depends  on  the  rate  of  shear  where  the  first  derivative  is  evaluated.  In 
general  protlems  the  speed  of  waves  depends  on  the  motion.  The  second 
order  general  problems  the  speed  of  waves  depends  on  the  motion.  The 
second  order  modulus  is  an  instantaneous  evaluation  the  second  Frechet 
derivative  at  the  wave. 
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The  amplitude  equation  of  Coleman  and  Curtin  is  notable  because  it 
is  simple,  general,  rigorous  and  Implies  interesting  physical  results. 
They  showed  that  an  initial  Jump  discontinuity  may  either  decay  or  grow 
depending  on  the  sign  and  magnitude  of  the  initial  discontinuity.  The 
assumed  Jump  discontinuity  may  blow  up  in  finite  time.  The  cause  of  this 
blow  up  is  associated  with  the  nonlinearity.  They  showed  that  an  accel¬ 
eration  wave  entering  into  a  region  at  rest  would  always  decay.  Thi3  is 
a  type  of  nonlinear  result  of  category  (il)  which  shows  that  nonlinear 
terms  can  lead  to  blow  up.  It  is  generally  assumed  that  the  loss  of  Cn 
regularity  (blow  up)  implies  the  formation  of  a  jump  in  the  n-1  deriva¬ 
tive  of  velocity  (shock  up).  The  decay  of  acceleration  waves  is  a  result 
of  category  (iii)  which  shows  that  nonlinear  terms  can  force  the  decay 
of  initial  discontinuities. 

It  was  not  clear  how  discontinuities  in  acceleration,  which  are 
equivalent  to  shocks  of  the  vorticity,  would  appear  in  the  fluid.  The 
results  discussed  in  the  appendix  by  M.  Slemrod  help  to  clarify  this 
issue. 

Slemrod  (1978)  showed  that  the  equation  of  motion  for  the  shearing 
perturbation  v(x,t)  of  a  shearing  motion  will  admit  a  differentiable  (in 
x.t)  solution  for  only  a  finite  time  for  appropriate  smooth  Initial  ve¬ 
locity  histories  when  the  constitutive  relation  is 


T<xy> 


e  t/A  vx(x,t-T)dt). 


(13.D 


/here  o  is  a  nonlinear  odd  function.  The  loss  of  differentiability  is 
issumed  to  imply  the  appearance  of  a  discontinuity  in  v,  as  shock  up  of 
,he  velocity,  a  vortex  sheet.  This  result  is  like  the  one  proved  by 
loleman  and  Gurtin,  but  implies  the  shock  of  smooth  data. 

In  the  appendix  of  this  paper  Slemrod  shows  that  when  the  constitu- 
ive  relation  is 


T<xy> 


eT/*  o(v  (x,t-t)dT 

0 


(13-2) 


or  a  nonlinear  o(’),  the  second  derivatives  of  v(x,t)  can  blow-up  when 
he  smooth  data  is  given  in  a  certain  way.  The  shock  up  assumption 
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which  one  is  obliged  to  associate  with  this  blow-up  of  second  derivatives 
is  the  appearance  of  a  Jump  in  vx,  a  vortex  shock  or  an  acceleration  dis¬ 
continuity.  It  is  interesting  to  note  that  T<xy>  given  by  0  3-2)  satis¬ 
fies  a  rate  equation 


±  T<xy>  ♦  iT<xy> 

dt  X 


o(vx(x,t)) 


(13.3) 


vaguely  resembling  some  popular  models. 

The  results  just  reviewed  show  the  remarkable  effects  of  the  choice 
of  constitutive  relations. 

It  is  possible  to  entertain  the  notion  of  successive  shock  ups  from 
smooth  data.  First  we  get  a  vortex  shock  from  nonlinear  effects  associ¬ 
ated  with  (13.2).  This  gives  an  acceleration  discontinuity  which  will 
lead  to  blow  up  of  the  vorticity  if  the  amplitude  of  the  vortex  shock  is 
larger  than  the  critical  amplitude  of  Coleman  &  Curtin.  In  reservation, 
I  wish  to  note  that  recent  calculations  with  popular  models  show  that 
the  critical  amplitude  can  be  Infinite. 

I  think  it  would  be  interesting  to  see  what  sort  of  blow  up  results 
could  be  obtained  for  shear  flows  of  Oldroyd  models  with  instantaneous 
elasticity. 

The  shear  flow  of  the  Oldroyd  models  (6.2)  are  governed  by  the 
system  of  first  order  quasilinear  equations. 


\  -  *  f  ■  °> 

't  *  [°-a)0  -  l]vx  *  x  ■  °- 


pv 


t 


(13.*0 


where  v(x,t)  is  the  rectilinear  velocity  in  the  direction  y  and 
(o,Y,t)  -  (T<yy>,  T<xx>,  T<xy>).  This  can  be  written  as 
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*  at  4  5  ax  ♦  2  ■  0 


he  characteristic  directions  are  given  by 

det|Ax  -  Bx  |  -  [p2xl  ♦  (1-a)o  -  n/l]  ■  0 

here  x  is  the  wave  speed.  The  streamlines  a^e  characteristic;  x  -  0, 
acause  streamlines  don't  move.  The  two  waves 

0 

x  - 


in  be  identified  with  waves  of  vorticity.  It  is  clear  that  the  speed  of 
le  wave  depends  on  how  the  liquid  is  stressed. 

We  should  like  to  know  if  smooth  initial  data,  given  to  (13.1!)  can 
Low  up.  This  depends  on  the  constitutive  equation,  as  we  have  already 
jen.  For  example,  the  problem  for  the  upper  and  lower  convected 
axwell  models  (a*-1)  cannot  lead  to  blow  up  of  smooth  data  because  the 
jverning  problems  are  linear. 

Johnson  and  Segalman  (1977)  have  shown  that  ( 1 3- •M) »  and  (13.1*)*  can 
‘  integrated  once: 
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>e  approximated  by  a  finite  series  of  multiple  integrals: 


t  (At(s)) 
3-0 


I 

n  odd 


Kp  (3i,..,sn) 


At(s,)...  At(sn)ds1...  dsn 


(2.3) 


restriction  to  odd  order  follows  from  the  isotropy  condition  (1.3).) 

assume  more  is  true,  namely  that  in  a  manner  similar  to  the  Taylor 
3  expansion  for  an  analytic  function  we  take  (2.3)  to  be  an  equality 
the  sum  may  be  infinite.  Choose 


(3| , .  •  • ,  3^) 


n 

-a  o  e 
n 


-a(s,+. 


.♦sn) 


1,3,5... 


a  is  a  positive  constant.  (This  is  consistent  with  fluid  of  fading 
y  type  [3],  [6].)  Finally  if  we  define 


*U)  -  £  0  £n 

* —  n 

n  odd 

itegrate  (2.3)  by  parts  we  obtain  (2.1). 

Constitutive  equation  (2.2)  is  equivalent  to  the  choice 


t  (At(s))  -  -  e-a3o(( At(s))  )ds  . 

s-0  Jo  3 

ternative  motivation  arises  from  consideration  of  a  material  of 
ype  where  Txy(t)  satisfies  the  ordinary  differential  equation 

£  Txy(t)  ♦  c«Txy(t)  -  o(vx(x,t))  (2.H) 


1  t,  -•  <  t  <  •  ,  and  all  x;  0  <  x  <  h  ,  where  a  as  before  is  a 
#e  constant.  Integration  of  (2.4)  in  the  usual  manner  for  first 
ordinary  differential  equations  yields  (2.2). 
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•  m 

t  (-A^(s))  -  -  t  ( Afc  (s))  , 

3-0  3*0 

m  m 

*i  (~At  (s))  -  ®i  (A*  (3))  ,  (i  -  1,2).  (1.3) 

3-0  3-0 

Here  subscripts  x,t  denote  partial  derivatives  with  respect  to  x  and  t 
respectively. 

A2.  CONSTITUTIVE  ASSUMPTIONS 

To  proceed  further  it  is  necessary  to  make  some  mathematical 
assumptions  as  to  the  nature  of  the  functional  t.  For  this  analysis  we 
assume  t  has  two  particularly  simple  forms,  i.e.. 


t  (A.ts)) 

s-0  1 


e  aavx(x,t-s)ds  , 


and 


(2.1) 


t  (A  (s))  -  e  aso(v  (x,t-s))ds  ,  (2.2) 

s-0  c  Jo 

where  a  is  a  real  valued,  odd,  analytic  function  defined  on  the  real  line 
and  a  is  a  positive  constant.  In  (2.1)  we  see  t  is  a  (generally)  non¬ 
linear  function  of  the  linear  functional 


I  e  asv  (x.t-s)ds  . 

Jo  x 

In  (2.2)  we  see  that  t  is  a  linear  functional  of  the  nonlinear  func¬ 
tion  a(vx(x,t-s)). 

A  motivation  for  the  choice  (2.1)  may  be  based  on  the  multiple  in¬ 
tegral  expansions  for  t  originally  presented  by  Green  &  Rivlin  [  H]  and 
Chacon  A  Rivlin  [5]  for  general  viscoelastic  materials.  In  [5]  the 
authors  showed  that  for  t  continuous  on  an  appropriate  function  space,  t 


>lds  with  Cl  replaced  by  C*.  D.  D.  Joseph  observed  that  this  may  have 
iportant  physical  consequences  as  now  the  appearance  of  Jump  dlscontln- 
ty  in  vx  may  be  expected,  i.e.  formation  of  a  vortex  shock. 

At  the  request  of  D.  D.  Joseph  I  have  prepared  this  short  note  ana- 
'zlng  the  different  effects  of  (1.1)  and  (1.2).  The  technical  details  of 
ie  breakdown  proof  are  omitted  and  may  be  found  [1]. 

.  RECTILINEAR  SHEARING  PLOWS 

If  in  a  fixed  Cartesian  co-ordinate  system  x,y,z,  the  velocity 
elds  of  a  flowing  fluid  body  has  the  form 

vx  -  0,  vy  -  v(x,t),  vz  -  0  (1.1) 

say  that  the  motion  is  a  rectilinear  shearing  flow.  For  such  a  flow 
e  condition  of  Incompressibility  div  v  -  0  Is  automatically  satisfied, 
leman  &  Noll  [3]  have  shown  that  if  the  fluid  is  a  simple  fluid,  then 
e  components  of  stress  obey  the  relations 

TXy(t)  -  t  ( At(s))  , 

3*0 

m 

TXX(t)  _  Tzz(t)  „  3l  (At(a))  f 
3-0 

CD 

Tyy(t)  -  TZZ(t)  -  3*  (Afc(s))  , 

3-0 

TXZ  -  TyZ  -  0  ,  (1.2) 

re  A*"  is  the  relative  shearing  history  defined  by 

•t 

At(s)  -  -  vx(x,t)dt  (0  i  s  S  «)  , 

•  t-s 

t.3,,3,  are  real  valued  functionals  obeying  the  identities 
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APPENDIX  A:  BREAKDOWN  OF  SHOOTH  SHEARING  FLOW  IN  VISCOELASTIC 
FLUIDS  FOR  TWO  CONSTITUTIVE  RELATIONS:  THE  VORTEX  SHEET  VS.  THE 
VORTEX  SHOCK. 

AO.  INTRODUCTION 

The  purpose  of  this  note  is  to  study  the  effect  of  choice  of  con¬ 
stitutive  relations  on  shearing  perturbations  of  steady  shearing  flows  in 
a  non-linear,  isotropic,  incompressible,  viscoelastic  fluid.  In  an  ear¬ 
lier  paper  [1],  I  showed  that  the  simple  constitutive  relation  relating 
the  shearing  stress  to  the  shear  rate  vx 


TXy(t) 


e  aTvx(x,t--r)dt). 


(0.1) 


yields  the  following  result: 

The  equation  of  motion  for  the  shearing  perturbation  v(x,t)  of  a 
non-trivial  steady  shearing  motion  will  admit  a  C1  in  (x,t)  solution  for 
only  a  finite  time  for  appropriately  chosen  smooth  initial  velocity  his¬ 
tories.  Here  o  Is  a  non-linear,  odd,  real  analytic  function,  o'  >  0,  a  is 
a  positive  constant. 

The  proof  of  the  result  hinged  on  the  fact  that  equations  for  the 
perturbed  fluid  motion  v(x,t)  can  be  written  as  a  non-linear  hyperbolic 
conservation  law  with  linear  damping.  In  this  form  the  equations  of 
motion  are  amenable  to  study  via  the  use  of  Riemann  invariants  and  an 
argument  of  Lax  [2].  As  a  side  benefit  from  this  formulation  continua¬ 
tion  theory  for  hyperbolic  equations  shows  that  loss  of  Cl  regularity 
implies  | v  |  |v  |  *  •  in  finite  time.  ;  This  suggests  the  appearance  of 

L  X 

A 

a  jump  discontinuity  in  v  ,  i.e.  formation  of  a  vortex  sheet. 

Sometime  after  appearance  of  my  result  C.J.S.  Petrie  remarked  to  me 
that  a  similar  result  could  be  obtained  for  the  choice  of  constitutive 
relation 


TXy(t) 


e  aTo(v  (x,t-i))dt 

0 


(0.2) 


for  a, a  as  above.  Upon  working  out  Petrie's  idea  I  noticed  the  results 
are  similar  but  not  the  same.  Specifically  the  above  mentioned  result 


| 
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have  been  smoothed.  On  the  other  hand  the  points  near  the  bottom  plate 
which  could  show  the  development  of  discontinuous  u/x  for  large  time  and 
large  W.  Perhaps  it  is  these  points  at  which  we  may  see  the  development 
of  shocks  of  u/x  are  only  what  PT&T  call  "unavoidable  noise."  It  is 
possible  to  find  pictures,  say  t  -  0.5,  W  -  1 ,  in  which  one  can  imagine  a 
point  of  discontinuity  of  the  slope  of  u/x  which  can  be  identified  with  a 
shock  in  the  vorticity. 

I  think  that  the  results  shown  support  the  notion,  already  explored 
in  Sections  10  and  12,  that  the  Weissenberg  number  is  a  measure  of  the 
fluids  elasticity  and  in  hyperbolic  problems  with  damping,  the  damping  is 
large  when  the  Weissenberg  number  is  small  and  vice  versa.  The  effect 
of  the  waves  is  almost  immediately  damped  in  the  case  W  -  0.1  and  is 
extraordinarily  persistent  when  W  -  500. 

It  would  be  interesting  to  see  if  and  what  discontinuities  develop 
in  the  squeeze  flow  problem  when  all  of  the  prescribed  data  is  smooth 
at  t  -  0. 
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An  initial  boundary  value  problem  for  the  reduced  quasilinear 
system  requires  that  one  set  of  boundary  and  initial  conditions.  At  the 
boundaries  y  -  0  and  y  -  H(t)  we  have  adherence 

u  »  v  »  0  aty-0, 

v  »  0,  v  «  -V  at  y  *  H(t). 

To  specify  initial  conditions  it  is  enough  to  give  the  fields  u  and  v, 
that  is,  f(y,t)  at  t  -  0.  PT&T  set  f(y,0)  »  -y.  Hence 

(u,v)  -  (x,-y)  (13.11) 

1 1-0 

which  is  an  lrrotational  squeeze  flow  with  vanishing  stresses  and  an 
outflow  independent  of  y  at  each  x.  The  wave  speeds  (13.8)  are 

dy/dt  -  C  at  the  bottom  plate 
1 1-0 

and 


dy/dt  -  -1  -  C  at  the  top  plate. 

1 1-0 

The  velocity  v  -  -y  at  t  -  0  is  compatible  with  the  prescribed  boundary 
conditions  but  u  -  0  on  y  »  0  and  on  y  -  H(t)  is  incompatible  with  u  -  x. 
The  velocity  v  and  all  the  stresses  are  continuous  as  t  +  0  but  the  ve¬ 
locity  u  »  0  at  the  walls  is  incompatible  with  (13.9).  In  a  linear 
theory  we  could  expect  this  Initial  discontinuity  to  propagate  into  the 
Interior  with  a  decay  in  the  magnitude  of  the  jump  in  u.  This  does  not 
appear  to  be  true  in  the  results  presented  by  PT&T. 

In  Fig.  13.2  I  have  reproduced  their  results  for  the  evolution  of 
-fy(y,t)  -  u/x  for  R  -  1  and  increasing  W.  The  graphs  of  the  evolution 
of  f  -  v  do  not  show  Interesting  features  and  are  not  presented  here. 

I  cannot  find  evidence  for  propagating  discontinuities  of  u/x  near 
the  top  wall.  The  wave  in  u/x  propagating  from  the  top  plate  seems  to 


I 


I 


W  -  Vl/h,  la  the  Weissenberg  number,  B  la  diagonal  with  entrlea 
[0,0,W,W,W,W,R]  and  D  dependa  on  g  but  not  on  lta  derlvatlvea.  The  wave 
apeeda  for  thla  system  y  ■  dy/dt  are  given  as  the  real  roots  of 

det[A  -  yB]  -  W^f-yJtM^f-y)*  -  1  -  WY]  -  0 

where  M*  -  RW.  The  normal  component  of  velocity  v  -  f  is  triply  charac¬ 
teristic.  The  wave  speeds  for  the  vortlclty  are  given  by  the  dimension¬ 
less  form  of  (13.8)« 


y  -  f  ± 


Numerical  results  of  PT&T  show  that  Y  ♦  0  as  W  ♦  •  in  such  a  way  that  1 
♦  WY  >  0.  There  is  no  Hadamard  instability. 


w’  .  W  c.*.  -L.’.  > 
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2 1 
dt 


0  3.8) 


where  C  -  /n/Ap  is  the  speed  of  small  amplitude  vorticity  waves  into 
regions  of  uniform  motion.  The  wavespeed  formula  (13.8)  shows  how  the 
speed  of  vorticity  for  solutions  of  the  type  0  3.6)  depends  on  the  point- 
wise  values  of  the  normal  component  of  velocity  and  the  tensile  stress  Y 
-  iyy.  Every  problem  which  can  be  represented  in  03-6)  will  have  wave 
speeds  for  the  vorticity  which  will  satisfy  (13.8).  Numerical  results  of 
PT&T  suggest  that  for  their  problem  Y/p  ♦  C1  >  0. 

In  principle,  it  should  be  possible  to  determine  some  of  the  effects 
of  nonlinearity  on  propagation  by  more  careful  studies  of  the  problem 
treated  by  PT&T. 

Let  h«  be  the  initial  distance  between  the  plates  and  let  v  be  the 
constant  squeezing  velocity 

(y.t)  -  r-(y.vt) 
n0 

are  dimensionless  variables,  also  called  y  and  t.  The  configuration  of 
flow  is  shown  as  Fig.  13.1.  The  dimensionless  velocity  components  are 
given  by 


(u.v)  -  (xfy  (y,t),f(y,t) 


and  have  been  made  dimensionless  with  V.  The  stresses  are  made  dimen¬ 
sionless  with  nV/h0. 

The  dimensionless  stress  and  pressure  it  are  then  represented  by 


r  <xx>  <yy>  _<xy>  „  1  , 

Ct  .  t  ,T  ,P]  -  X,  ♦  x%,xT,Y,  2X  Pi  +  Pa 


03.9) 
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I  have  analyzed  their  problem  using  some  elementary  hyperbolic 
theory.  I  work  with  their  reduced  equations.  The  reduced  equations 
given  by  (13.10)  are  a  consequence  of  representation. 


[u,v]  -  V[-xfy(y,t),f(y,t)] 


(13.6) 


Here,  subscripts  denote  differentiation.  We  treat  the  rest  of  the  reduc¬ 
tion  later.  For  now,  consider  Eq.  (6.4)  with  a  -  1  for  the  vortlcity 


p  0  +  pu  &k  *  2pv  +  (pu2"°'a)If 


+  2(Puv+t)^  ♦  (PV*-T-J)0  -  L.O.T. 


(13.7) 


The  vorticity  compatible  with  (13.6)  satisfies 


r  .  -1  *  -  V  Xf 
4  2  3y  2  yy 


def 


XU, 


C  *  u,  c  «  0. 

^x  ’  ^xx 


These  relations  imply  that 


32c  1  [*3u  3uj1 

[3x3t*3x3yJ  *  [_3t ’  3y J 


are  first  derivatives  of  u  and  hence  of  lower  order.  We  may  then  write 
(13.7)  as 


A 


3t* 


♦  2B 


32C 

3y3t 


L.O.T. 


where  A  -  p,  B  -  pv,  C  -  pv2  -  Y  -  j.  The  characteristic  equation  for 
this  second  order  problem  are  given  by 
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<(s)  -  vx(x,s) 

is  the  shear  rate  at  the  layer  x  at  time  s  i  t.  It  follows  that 


I«t> 


't 

t 

1  -  /I -a5 

e-(t-s)/X  K(a)3in(/f^ 

K(s')ds')ds 

-• 

s 

We  may  consider  fluids  with  short  memories,  small  1.  After  Integrating 
by  parts  we  find  that 

tt  -  -J  x(t)[1  -  (1  -  a*)X*K*(t)  ♦  OCX*)]  ♦  j  -  0  0  3.5) 

This  is  in  the  form 

^  T<Xy>  +  I  T<Xy>  “  X  vx(x,t)[1  '  (1  *  a*)X*v*(x,t)  ♦  0(X,] 

which,  up  to  terms  0(X2)  is  like  03.3),  leading  to  the  blow  up  of  vorti- 
city  when  1  -  a2  <  0. 

We  turn  next  to  a  still  more  complicated  flow  involving  wave  propa¬ 
gation,  which  is  originally  set  in  two  space  dimensions. 

N.  Phan  Thien  and  R.  I.  Tanner  (1983).  hereafter  called  PT&T,  have 
considered  the  problem  of  flow  of  an  upper  convected  Maxwell  model 
which  is  induced  by  squeezing  the  fluid  between  infinite  parallel  planes. 
They  do  not  linearize  or  neglect  terms.  They  call  their  solution  exact 
because  they  find  a  separable  solution  03*6)  which  reduces  the  governing 
three  dimensional  qqasilinear  system  to  a  nonlinear  two-dimensional 
system  which  they  integrate  numerically.  They  have  impulsive  initial 
data  in  the  following  sense.  Initially  an  irrotational  plug  flow  is  pre¬ 
scribed.  The  velocity  v  normal  to  the  plates  increases  linearly  from 
zero  at  the  bottom  to  the  prescribed  value  at  the  top.  At  this  same  in¬ 
itial  Instant  the  velocity  u  parallel  to  the  plate  is  squeezed  out  as  a 
plug  flow.  Since  u  must  vanish  at  the  plates,  this  component  is  pre¬ 
scribed  as  discontinuous.  In  linear  problems  this  discontinuity  would 
propagate  and  decay  along  characteristics.  PT&T  do  note  that  their 
numerical  work  shows  wave  propagation  but  their  graphs  are  confusing;  it 
is  hard  to  tell  what  is  propagating  and  how  it  is  propagating. 
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A3.  SHEARING  PERTURBATION  OP  A  STEADY  SHEARING  PLOW 

Let  us  assume  a  viscoelastic  fluid  satisfies  either  the  constitutive 
relation  (2.1)  or  (2.2).  Consider  the  problem  when  the  fluid  is  confined 
between  two  parallel  walls  of  infinite  extent  at  x  -  0  and  x  •  h.  The 
top  wall  at  x  -  h  moves  with  velocity  V.  In  the  absence  of  body  and 
driving  forces  the  equation  of  onservation  of  linear  momentum  (see 
Coleman  &  Curtin  [7])  becomes 


pvt(x,t) 


e  03vx(x,t-s))ds)x 


for  constitutive  relation  (2.1)  and 


(3-D 


e  a3c(v  (x,t-s))  ds  (3.2) 

0 

for  constitutive  relation  (2.2).  Here  p  is  an  (assumed)  constant  mass 
density,  p  >  0  .  Also  we  consider  the  case  of  no-slip  boundary  condi¬ 
tions 


pvt(x,t)  - 


v(0,t)  -  0  ,  v(h,t)  -  V  .  (3-3) 

Systems  (3-D*  (3>3)  and  (3.2),  (3.3)  admit  the  steady  rectilinear 
flow  (vt(x,t)  -  0)  solution 

V 

vw  -  -a  . 

To  study  stability  of  the  flow  against  shearing  perturbations  we 
set 


v(x,t)  -  v(x,t)  -  • 

n 

We  observe  that  for  relation  (2.1)  this  implies 
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pvt(x,t) 


o( 


e_a3vx(x,t-3)d3  *  ±)x 


(3.“) 


with  boundary  condition 

v(0,t)  -  v(h,t)  -  0  (3.5) 

On  the  other  hand  for  relation  (2.2)  we  see  that  the  perturbation 
v(x,t)  satisfies 


pvt(x,t) 


e  a3o(vx(x,t-s)  ♦  ^)x  ds 

0 


(3.6) 


along  with  boundary  condition  (3-5). 

In  either  case  we  prescribe  a  smooth  velocity  history 

v(x,t)  *  v0(x,t)  -•  <  t  i  0 

consistent  with  either  (3.4),  (3-5)  or  (3-5),  (3.6).  Thus  the  two  fluid 
cases  are  governed  by  non-linear  boundary  initial  history  value  problems. 

A4.  ANALYSIS  OF  FLOW  WITH  FIRST  CONSTITUTIVE  ASSUMPTION:  STRESS 
NON-LINEAR  FUNCTION  OF  A  LINEAR  FUNCTIONAL  OF  SHEAR  RATE. 

We  consider  the  perturbed  flow  governed  by  constitutive  relations 
(2.1).  In  this  case  the  time  evolution  is  governed  by  (3.^)»  (3.5).  To 
simplify  matters  we  write 

k0  '  j[0('  *  3T>  '  "‘SS’]  • 


u(x,t) 


—  Q«  A 

e  v  (x,t-s)ds  , 

0 


—  Ml  A 

a  vx(x,t-s)ds  . 


w(x,t) 


we  see  o  satisfies 
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o(0)  ■  0  , 


o"(0)  *  0  . 


We  further  Impose  the  condition  that  (4.2)  be  strictly  hyperbolic, 
i.e.  the  matrix 


[:  ;1 


possesses  real  distinct  eigenvalues.  Strict  hyperbolicity  is  easily  seen 
to  be  equivalent  to  the  condition  o’  >  0  . 

A5.  ANALYSIS  OP  FLOW  WITH  SECOND  CONSTITUTIVE  ASSUMPTION:  STRESS 
LINEAR  FUNCTIONAL  OF  A  NON-LINEAR  FUNCTION  OF  SHEAR  RATE. 

We  consider  the  perturbed  flow  governed  by  constitutive  relation 
(2.2).  In  this  case  the  time  evolution  is  governed  by  (3-5),  (3.6).  Again 
to  simplify  matters  we  write 


.  V,  V 
*  h>  -  •  h  • 


u  ■  vt  • 


We  note  (3.6)  can  be  written  as 


t 

\ 

.  —  • 


vt  -  e  a(t-T)  o(vK(x,t»x  dT 


*  : 


which  upon  differentiating  with  respect  to  t  yields 


vtt  "  ~a  e  a(t  Z)  o(^x  (x*T))x  dT  +  a(vx(x,t))x 
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vtt  ♦  avfc  -  °(vx(x,t))x  (5.B) 

where  we  have  used  (5.1). 

From  the  definitions  of  u,  w  above  we  see 

"t  "  “x  * 

“t  “  o(w)x  -  ou  ,  (5.4) 

with  boundary  conditions 

u(0,t)  -  u(h,t)  -  0  ,  (5.5) 

and  initial  conditions 
u(x,0)  -  u0(x), 

w(x,0)  -  w0(x),  0  <  x  <  h  . 

Again  the  values  of  u0(x),  w0(x)  can  be  obtained  from  the  respective 
definition  of  u,w  in  terms  of  the  given  velocity  history  v0(x,t), 
*•  <  t  S  0. 

Also  we  note  that  if  o"(0  *  0  for  some  c  real  and  when  the  speed 
of  top  wall  is  given  by 

V  -  he 


then 


o(0)  -  0  , 

o"(0)  *  0  .  (5.6) 


As  before  the  condition  o'  >  0  implies  (5.4)  is  strictly  hyperbolic 
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A6.  A  BREAKDOWN  RESULT. 

We  have  shown  in  sections  ^  and  5  that  evolution  of  perturbed  flow 
is  governed  by  a  system  of  the  form 


W 


t 


U 


x 


Ut  ■  K(w)x  '  aU  * 


(6.1) 


with 


U(0,t)  -  U(h,t)  -  0, 

l(x,0)  -  U0(x),  W(x,0)  -  W0(x)  (6.2) 


and 


K’  >  0,  k(0)  -  0,  K"(0)  4  0. 

Analysis  of  (6.1),  (6.2)  has  been  given  in  [1].  We  shall  not  repeat  that 
analysis  but  only  state  the  relevant  breakdown  result. 

Define  Riemann  invariants  for  (6.1)  by 

r  -  U  ±  0(W)  (6.3) 

s 

where 

fw  _ 

0(W)  -  /K'(s)  ds. 

Jo 

The  transformation  given  by  (6.3)  from  (U,W)  R  x  R  to  (r,3)  R  x  R  is 
one-one.  Also  we  assume  the  initial  data  r(0,x)  ■  r0(x),  s(0,x)  -  s0(x)  to 
be  smooth  functions. 

Our  main  breakdown  result  is  as  follows. 

Theorem  6.1:  Suppose  |r0|,  |s0l  are  sufficiently  small  and 

K'(0)  >  0,  K"(0)  >  0.  If  r0,x  or  s0,x  is  positive  and  sufficiently  large 
at  any  point  x,  then  (6.1),  (6.2)  has  a  solution  (W,U)  in  C‘[0,h]  x  C'[0,h] 
for  only  a  finite  time.  A  similar  result  holds  if  K"(0)  <  0  and  r0,x  or 


s0,x  is  sufficiently  negative  at  any  point  x. 

Proof:  See  the  proof  of  Thm.  3*1  in  [1]. 

Standard  existence  theorems  (see  for  example  Section  1.8  of  [8]  and 
Chaps.  2  and  3  of  [93)  imply  that  under  the  hypothesis  of  Theorem  6.1  we 
have  |Ut|  ♦  |Ux|  ♦  jwj  ♦  |Wx|  ♦  •  in  finite  time.  This  suggests  but 
does  not  prove  the  occurrence  of  a  jump  discontinuity  in  U  and  W,  l.e. 
the  formation  of  a  shock. 

Also  we  note  that  since 

Ix  -  Ux  *  ♦'<«"* 

we  will  have  r0j{  large  if  or  is  large,  s0j(  large  if  or  -W^  is 
large,  -r„x  large  if  -% or  -W^  large,  - s ^  large  if  -U^  or  W^  large. 

7.  PHYSICAL  IMPLICATIONS  OP  BREAKDOWN  OP  SMOOTH  SOLUTIONS 

We  examine  the  implications  of  the  breakdown  result  of  Section  6 
with  respect  to  our  two  constitutive  relations. 

First  we  consider  constitutive  relation  (2.1)  In  this  case  we  see 
that  if  o"(V/ah)  ^  0  and  u,x(x)  or  w#x(x)  is  appropriately  sufficiently 
large  of  sufficiently  negative  u^  (depending  on  the  sign  of  o"(V/ah)), 
Theorem  6.1  Implies  |ux|  ♦  |wx|  ♦  »  in  finite  time. 

Since  we  know  from  (3.6)  that 


and  from  the  definition  of  u  we  have 

\  “  ux  *  <**• 

| u  |  ♦  | w  |  ♦  •  in  finite  time  implies  Iv  I  ♦  I v  |  ♦  •  in  finite  time. 
Again  this  suggests  but  doesn't  prove  that  v  and  hence  v  forms  a  Jump 
discontinuity  in  finite  time.  In  this  case  the  singular  surface  across 
which  the  discontinuity  in  v  is  called  a  vortex  sheet. 

Next  we  consider  the  case  of  constitutive  relation  (2.2).  In  this 
case  we  see  that  if  u«x(x)  •  vtx(*t0)  or  w#J((x)  -  vxx(x,0)  is  appropri¬ 
ately  sufficiently  large  or  sufficiently  negative  (depending  on  the  sign 
of  o"(V/h))  Theorem  6.12  implies  | vfcx |  ♦  |vxx|  ♦  •  in  finite  time.  Once 


more  this  suggests  but  doesn't  prove  that  either  v.  or  v  forms  a  Jump 

v  X 

discontinuity  in  finite  time.  However  from  the  Rankine-Hugoniot  Jump 
condition  for  (5. 4) 


-(ff)Cw]  -  Cu] 

-(^[U]  -  Co(w)]  (7.1) 

at 


we  see 


-(3?*V  -  K  -  K>]  ■ 

where  x  -  s(t)  denotes  the  surface  across  which  the  Jump  occurs.  Hence 
if  [v  ]  4  0  then  [o(v  )]  4  0  (since  o'  >  0)  and  by  (7.2)  [v  ]  4  0.  Con- 

XX  v 

versely  if  [vt]  4  0  then  from  (7.1)  we  have 

-<£>[.„:  -  [vt]  (7.3) 


and  hence  [v  ]  4  0.  Thus  the  appearance  of  a  jump  in  v  or  v  implies  a 

X  v  X 

Jump  in  the  other. 

We  define  a  propagating  singular  surface  across  which  the  accelera¬ 
tion  vfc  experiences  a  jump  discontinuity  as  an  acceleration  wave.  Simi¬ 
larly  we  define  a  propagating  singular  surface  across  which  the  vorti- 
z  y  z 

city  w  -  curl(v  ,z  ,y  )  »  v  (x,t)e  experiences  a  jump  discontinuity  as  a 
vortex  shock.  Our  analysis  shows  that  for  our  flow  the  vortex  shock  and 
acceleration  waves  are  equivalent  and  can  be  expected  to  form  in  finite 
time  if  constitutive  relation  (2.2)  holds. 

We  thus  see  the  remarkable  effect  of  choice  of  constitutive  rela¬ 
tion.  In  one  case  (2.1)  appropriately  chosen  initial  data  appears  to 
force  the  formation  of  a  Jump  in  v  (a  vortex  sheet)  while  in  the  second 
case  (2.2)  appropriately  chosen  initial  data  suggests  formation  of  a  Jump 
discontinuity  in  v  (a  vortex  shock  and  equivalently  in  v  (an  accelera- 
tlon  wave). 
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